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The absolute intensities of the chief infra-red absorption bands of ethylene and nitrous oxide 
have been determined experimentally to be (in cycles per centimeter at N.T.P.): 


Ethylene 
950 cm™ band: 1540 10" 
1444 cm band: 21010" 
3050 cm! band: 84010! 


Nitrous oxide 
590 cm! band: 12010" 
1285 cm band: 115010" 
2224 cm™ band: 5600 10" 


From these intensities it is possible to calculate the dipole moments, yu, of bonds between 
vibrating atoms and also the rates of change of dipole moment with inter-nuclear distances, 
Ou/dr. The values of uw obtained are all of the order of magnitude of 0.5107" c.g.s. unit, 
while values of 0u/0r are of the order of 1X 10-" c.g.s. unit. The results are in agreement with 
available data on infra-red dispersion and on atomic polarization within the limits of experi- 
mental error, which are estimated to be approximately +10 percent. 





N a recent paper! a method has been described 
for measuring the intensities of infra-red 
vibrational absorption bands. The integrated 
apparent absorption coefficient per molecule is 
measured at several partial pressures of absorbing 
gas and extrapolated to zero partial pressure. 
This limit can normally be expected to equal the 
true integrated absorption coefficient. The use 
of a non-absorbing gas to broaden the individual 
lines composing the band assures that the limit 
will be approached closely for measurable values 
of the absorption. This pressure broadening need 
not be so great as to merge completely the 
rotational lines, but must reduce the variations 
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of a(v) across the spectrometer slit until the 
observed absorption by a given amount of gas 
is not appreciably increased by a further increase 
of non-absorbing foreign gas pressure. It has 
long been known that the addition of a non- 
absorbing gas causes, in some cases, an increase 
in the measured absorption. The increase caused 
by the first increment of foreign gas is consider- 
able, but further additions produce smaller and 
smaller effects.2 The pressure required in practice 
varies from zero to about 50 atmos., depending 
primarily upon the spacing between the rota- 
tional lines of the absorbing gas. 

The experimental procedure is simplest for 
molecules whose lines are adequately broadened 
by one atmosphere total pressure. N2O and 


2 See, for example, E. von Bahr, Ann. d. Physik 29, 780 
(1909) ; 33, 585 (1910). 
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Fic. 1. Steps in calculating the intensity from experimental data (950 cm band of ethylene), 


CH, were chosen as fulfilling this requirement? 
without having unduly complicated spectra. This 
paper describes measurements of the intensities 
of their fundamental bands. 


EXPERIMENTAL PROCEDURE 


The absorption of the gas under study is 
measured as a function of frequency for various 
partial pressures of absorbing gas, at a fixed 
total pressure, and the intensity is obtained from 
the extrapolation previously described.! An auto- 
matic recording prism spectrograph was em- 
ployed® whose effective slit width varied from 
about 5 cm= to 25 cm depending upon fre- 

3 The instrument is that described by Gershinowitz and 
Wilson, J. Chem. Phys. 6, 197 (1938). Improved slits and a 
more sensitive thermocouple were used, and an off-axis 


parabolic mirror has replaced the original spherical one for 
the more recent of the measurements here discussed. 


quency. The sample was enclosed in a cell 30 cm 
long, with NaCl or KBr windows. Air was added 
as a broadening gas to a total pressure of about 
one atmosphere. The partial pressure of the 
sample was varied from nearly an atmosphere 
down to about 0.1 mm Hg for an intense band. 
The data for low pressures (that is, those which 
lead to a maximum observed absorption not 
greater than about 30 percent) are most im- 
portant in determining the intensity. The lower 
limit to the partial pressures employed was set 
by the fact that absorptions of less than about 
10 percent were too small to be determined 
precisely from the record obtained. 

The first measurements, made in 1940-41, 
have been repeated in 1946 with special emphasis 
on the low partial pressures. 

The basic data were obtained as a record on 
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photographic paper of galvanometer deflections 
for sample (7) and blank cells (JT) at frequency 
intervals averaging about 5 cm. The ratio of 
these deflections was then plotted as a function 
of frequency on semi-logarithmic graph paper, 
and the area between the curve for zero pressure 
of absorbing gas and that for pressure p was 
measured to give the integrated apparent absorp- 
tion, defined by (1), below: 


@=pLB= [ In(To/T)dv (1) 


where L is the cell length and » the frequency. 
The integration is over the band. 

While it is the integrated apparent absorption 
coefficient, B, that is extrapolated! to obtain the 
true intensity, A, this quantity is liable to 
fluctuate considerably if calculated directly from 
experimental data. It was most convenient to 
plot the integrated apparent absorption, ®, draw 
a smooth curve through the experimental points 
and select smoothed values of @ from which to 
calculate B for extrapolation. These steps are 
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illustrated by data for the 950 cm band of 
CoH, in Fig. 1. 


EXPERIMENTAL DATA 


The integrated apparent absorption, ®, is ob- 
tained from the basic data by unambiguous 
arithmetical reduction, whereas quantities ob- 
tained through smoothing and extrapolation 
involve individual judgment. The plotted values 
of @ in Figs. 2 and 3 are, therefore, a summary 
of the basic data, but the values of the intensities 
given in Table I involve judgment as well. 

Data for NO are presented in Fig. 2, and 
those for CeH, in Fig. 3. The case of NO is 
clear-cut, but that of CsH, has some unsatis- 
factory aspects. In the first place, 1941 and 1946 
data are not in agreement for the 950 cm band, 
and the discrepancy is more than the expected 
limits of experimental error (see discussion 
below). In addition, both 3050 cm and 950 cm 
“bands” actually involve absorption by two 
vibrational modes which are not resolved. The 
separate contributions can only be estimated in 
these cases. 
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Fic. 3. Experimental results for CoH. 


ACCURACY OF RESULTS AND LIMITATIONS 
OF METHOD 


The most serious general limitation of the 
method is the critical importance of the low 
pressure measurements where a small absolute 
error results in a large percentage error. In- 
accuracies in measuring records, and in deter- 
mining the ratio of sample to blank deflections in 
the absence of absorbing gas, may be responsible 
for errors of 5-10 percent. By taking a very 
large number of measurements the error in the 
smoothed curve could be made less than this, 
but it is believed that the intensities given in 
Table I may be in error by about 10 percent. 

A number of other sources of error may be 
minimized by proper procedure. Stray light must 
be corrected for, and this correction becomes a 
serious potential source of error if the contribu- 
tion of stray light to the deflections is 20 percent 
or more. In the above data only the N2O 590 cm= 
band was low enough in frequency to have 
appreciable stray (about 40 percent). Rapid 
fluctuations in the intensity of the incident beam, 
due to absorption by atmospheric CO2 and H,0O, 


are a potential source of error,! even if accurate 
ratios of transmission through the filled and 
blank cells are obtained. During the 1946 meas- 
urements an attempt to eliminate this absorption 
was made in the case of the 1444 cm band of 
C2H,, the only one of the bands studied which 
lies in a region of strong atmospheric absorption. 
While the H.O absorption was by no means 
eliminated, it was much reduced. Data on this 
band show good agreement between 1941 and 
1946 data, suggesting that the errors due to 
atmospheric absorption are not serious in this 
case. 

A final source of error which is difficult to 
eliminate is that of ‘‘envelope effect’ in the 
Q-branches of perpendicular bands, since these 
branches may be intense and narrower than the 
slit width of low resolving power spectrographs. 
To give the Q-branch a half-width of 10 cm by 
pressure broadening (assuming it to be effectively 
a single line) would require about 100 atmos- 

*One would expect atmospheric absorption to have no 
effect if randomly distributed in frequency, but to have a 


serious effect if the lines are spaced regularly in relation to 
those of the absorbing gas. 
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pheres broadening gas pressure. To eliminate 
the envelope effect by extrapolation to zero 
equivalent path length may require measure- 
ments at impracticably low absorptions. This 
type of difficulty may well account for the 
difference between 1941 and 1946 data on the 
950 cm= band of CsHg,, since this band is the 
only one of those measured which has prominent 
multiple Q-branches (spaced about 10 cm 
apart), and the resolution of the spectrograph 
was considerably better in 1946 than in 1941. 

In order to ascertain whether the resolution of 
the spectrograph in 1946 (about 10 cm effective 
slit width in this region) was adequate, a series 
of measurements was made on an instrument of 
considerably higher resolving power (a Perkin 
and Elmer spectrograph at Brown University).§ 
While the fine structure was very much more 
pronounced on these records, the values of 
integrated apparent absorption were in good 
agreement with our 1946 data, and the extrapo- 
lated intensities very nearly the same. We may 
conclude that errors due to envelope effect are 
not serious in this case. 

On the basis of these estimates of errors, it is 
considered that the results above may be in 
error by about 10 percent, but that certainly in 
no case should errors exceed 25 percent. In order 
to achieve this accuracy, considerable care must 
be employed in measuring ‘the records and a 
sufficient number of pressures must be used. 

The preliminary results,® published in 1941, 
were satisfactory for N2O but were based on too 
few points to give good agreement with our 
final results on ethylene, even when adjusted 
to N.T.P. 


INTERPRETATION OF RESULTS 
The integrated intensity A; is related to the 
matrix elements of the electric moment by 


A;= | a(v)dv=(82*Nv;/3hc) 


xX { | win | 2+ | prey’ 9] 2+ | pecs | 2}, (2) 


in which the integration covers the ith band, of 





° We are indebted to Professor Paul C. Cross, Professor 
Donald H. Hornig, and Mr. Earl Wilson of Brown Uni- 
versity for their cooperation with these measurements. 
(1944) J. Wells and E. B. Wilson, Jr., J. Chem. Phys. 9, 659 





frequency »;. a(v) is the true absorption coeffi- 
cient, NV is the number of molecules per cc, h is 
Planck’s constant, and c the velocity of light. 
Miz” etc., are the matrix elements of the com- 
ponents of the electric moment for a transition 
between the ground state and first-excited state 
of the ith normal mode.’ This quantity may be 
calculated on the basis of various simple models. 

If, for example, we assume the wave equation 
to have been separated in normal coordinates, 
Q;, we can, to a first approximation, express the 
(vector) dipole moment, u, as a linear function 
of the 0;: 


u—wo= Di (dn/0Q,)0i, (3) 
ui®°= (du/0Q;)Q;"°. (4) 


Q;° is the matrix element of Q; for the 0-1 transi- 
tion in question; wo is the constant part of the 
dipole moment. Each normal coordinate may be 
treated as a harmonic oscillator coordinate for 
which 


whence 


Qi*°= (h/8x?»4)}. (5) 


Substituting (4) and (5) in (2) yields a simple 
relation between the observed intensity and the 
coefficient in Eq. (3): 


A;=(Nr/3c) { | Ouz/dQ;|* 
+ | Ouy/9Q;|?+ | du./dQ;|?}. (6) 


™ In order to explain the observed values of the 
A; with a smaller number of constants we might 


TABLE I. Intensities obtained by extrapolation. 











Intensity 
(cycles per centimeter 
Substance Band at N.T.P.) 
NO 590 cm7 120 10” 
NO 1285 cm™ 1150 10" 
N.O 2224 cm 5600 x 10° 
CH, 950 cm7 1540 X 101% b 
CH, 1444 cm™ 210 10" 
CoH, 1890 cm! (combination) 67 x 10" 
CH, 3050 cm7 840 x 101% 








® The intensity of 1540 X10" for the 950 cm=! “‘band” is based on 1946 
data only, On account of better resolving power it is considered to be 
more accurate than that obtained in 1940-41. 

b This includes a contribution from the weak 995 cm~ band. From the 
observed asymmetry of this pair of overlapping bands we estimate that 
the total intensity should be apportioned as follows. 950 cm=}: 1390 
X10 cycles; 995 cm=!: 150 X10" cycles. ‘ 

© There are two overlapping bands here, at 2990 and 3105 cm~1, Using 
the relative intensities observed by Levin and Meyer, J. Opt. Soc. Am. 
16, 137 (1928), we tentatively divide this observed intensity as follows. 
2990 cm~!: 270 X10" cycles; 3105 cm~=: 570 X10" cycles, 


7G. Herzberg, Infrared and Raman Spectra of Polyatomic 
— (D. Van Nostrand Company, New York, 1945), 
p. 261. 
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TABLE II. Normal coordinates for ethylene. 











Sym- Fre- 
Mode metry quency* Normal coordinate 
vn Bs, 2989.5 Qu = 12.6 10713.S,—1.95 « 1076S; 
Vi2 Bou 1443.5 Qe =—1.12X 10-,5,+-90.8 X 10-6S; 
vy Bo, 3105.5 Qg=—12.2K10-%S,+-0.72 xX 10-6S, 
v0 §6BBay 995 Qio = +0.91 10735,+-93.1 x 10-6S, 
vr By 949.2 Q;=126X10-*Si 








* See reference 7. 


endeavor to reproduce the intensities by assign- 
ing fixed “effective charges’’ to each type of 
atom, and taking the form of the normal modes 
into account. Such a simplification would be too 
drastic, however, since the intensities calculated 
in this manner do not agree well with those 
observed. It is more satisfactory to assume, as 
has been done, for example, by Rollefson and 
Havens,* and Wolkenstein,® that the electric 
moments of the bonds are additive. Then the 
intensities of the stretching vibrations depend 
on du/dr, and those of the bending vibrations 
depend on uw, r being the bond length. 

This relation can be written explicitly since 
the normal coordinates, Q;, are linearly related 
to the internal coordinates, R, (e.g., changes in 
inter-atomic distances, etc.), by the transfor- 
mation 


Re= Di LiiQi, Qi= Le (ZL) aRe. (7) 


But the dipole moment of the molecule is related 
to the coordinates, Q;, by the du/dQ; as in (3). 
From (3) and (7) we obtain 


u—wo= Di (0u/90:)0: 
= Dix (0n/0Q;)(L) x Re 
= Di (Ou/AR:) Re, 
du/OR.= Di (Ou/0Q,) (L) x. (8) 


Since the transformation L~ is easy to normal- 
ize,’ this is a convenient method of determining 
the behavior of bond dipole moments with 
respect to stretching, bending, or other internal 
coordinates from the observed absorption in- 
tensities."! 


8 R. Rollefson and R. Havens, Phys. Rev. 57, 710 (1940). 

®M. W. Wolkenstein, J. Phys. USSR, 5, 185 (1941). 

1 FE. Bright Wilson, Jr., J. Chem. Phys. 9, 76 (1941). 

11 For symmetrical molecules it may be convenient to 
introduce internal symmetry coordinates, ®, as an inter- 
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TABLE III. Experimentally determined characteristics of 
bond moments of C2H, (c.g.s. units). 











Bond Bands involved im u/s Ou /dr 
CH v1: out-of-plane bending 0.77 X10—18 0.72 X10- — 
CH  v»v7in-plane bending 0.52 X107-8 0.48 X10-10 — 
CH _ »s in-plane bending 0.37 X10—18 0.33 X10- — 
+0.70 x 10-1 
CH vw stretching a = or 
—0.60 10-1 
+0.57 K 10-10 
CH ve stretching — —_ or 
—0.53 K10-10 








This procedure involves one serious difficulty, 
however, arising from the fact that it is the 
squares of duz/0Q;, Ou,/0Q;, and du,/dQ; that 
enter into (6). The signs of these coefficients are, 
therefore, undetermined, and can be either posi- 
tive or negative. When they are substituted in 
Eq. (8) a large number of combinations of signs 
are possible, each leading to a different value of 
du/dR;,. Other arguments must be introduced 
in order to eliminate as many of them as possible. 

Following this procedure bond moments and 
their rates of change were calculated from the 
intensities given in Table I. In order to obtain 
the (L~') coefficients normal coordinate treat- 
ments of N2xO and CH, were carried out using 
the methods previously described.” 


Ethylene 


In the case of €2H, the normal coordinates 
were calculated from the observed frequencies on 
the assumption that the dependence of the 
potential energy upon the infra-red active sym- 
metry coordinates can be expressed without 
interaction terms as 


V=3 i fuSetfhosSetfeoSetfraS7?t+fu S17}, 


TABLE IV. Normal coordinates for N2O. 











fi fis fo (L)1 (L~!)12 (L-)or = (LZ )at 
in units of 105 in units of 10-12 
0:50 15.47 14.05 —2.30 +165 +3.24 +3.78 
1.36 17.91 11.36 —2.72 +1.13 +2.89 +3.95 
2.72 20.32 10.29 —3.03 +0.613 +248 +4.07 
5.00 22.94 9.884 —3.42 -—0.017 +2.05 +4.22 
10.00 26.51 11.38 —3.85 —1.14 +1.05 +4.01 








® There are, of course, two algebraic solutions for fi and f2 with each 
value of fiz assumed, of which only one is tabulated. The other leads, in 
each case, to a solution with fi smaller than f2, and was discarded since 
the 2}-fold NN bond must surely be stiffer than the 14-fold NO bond. 


mediate step, determining first the du/@® as in (8) and 
then the du/dR;. 
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TABLE V. Characteristics of bond moments of NO (c.g.s. units). 

















fie (dn/8r NN (du/8r) NO -NN ¥NO UNN/S #NO/S 
1.0 —2.9X 10-1 —7.5X10-% 
2.0 —4.110-" —6.6X 10-1 
= 10 am —10 
oe aint «= ~aeeip-op 050X107" = —0.67X10-*# = -0.45x10-" = —0,56x 107° 
7.0 —7.1X 107% —3.4X 10-1 
10.0 — 8.2 10- — 1.6 10-1 
where opposite signs, leaving two possibilities which 
Keticinmte~ee lead to the two values of du/dr for vy, and v. If 
er P : }aton C—C axis the 0u/0Q’s had the same signs, these two CH- 
Ss=3rRiyitrve—-vs—73}, stretching vibrations would involve values for 
du/dr of opposite sign, which seems quite im- 
a? _— ’ 
Ss=F(ri—retrs—ra}, | Perpendicular to possible. 
S7=3rR\ yi-yv2t 73-74}; C—C axis 


Siu=(L(—rR cosy) /v2]{5a+ 6}. 


The internal coordinates r; and y; are changes in 
CH bond lengths and C—C—H bond angles, 
respectively, and 6, and & are out-of-plane 
bending angles of the H—C—H planes with 
respect to the C—C bond. 

The force constants were determined to be 


fas=5.07 X 10° dynes/cm, 

fss=6.33 X 10” dynes/cm*/radian?, 

fos =5.08 X 10° dynes/cm, 

fiz=3.07 X 10” dynes/cm*/radian?, 
fir 1 =5.11 X 10” dynes/cm*/radian?. 


The corresponding normal coordinates are 
given in Table II. 

The results for CH, are presented in Table III. 
Values are given for u, 0u/dr, and u/s, where s is 
the equilibrium inter-atomic distance. (It is 
assumed throughout that the dipole moments are 
directed along the bond in question.) u/s and 
du/dr have the dimensions of charge, and may 
be compared with the charge of —4.8X10~° for 
an electron. 

In Table III considerations of the electro- 
negativities of the atoms have been employed to 
determine the signs when possible. Since hydro- 
gen is less electronegative than carbon, negative 
values for the CH bond moment would be 
unreasonable, and this restriction allows us to 
select the proper sign for three of the five 
0u/0Q’s for C2Hy. The remaining two must have 


Nitrous Oxide 


The case of N2O is made more complicated by 
the existence of resonance between the two 
structures 


(a) N=N=O, (b) N=N-O. 
The potential energy must be expressed as 
V=${firet2friretfor?t+fsy’}, 


where 7;=NN stretch, 72=NO stretch, and 
y¥=NNO bending. The interaction coefficient 
cannot be evaluated from the two observed 
frequencies, but is probably large (and positive) 
since a stretch of the NO bond will favor struc- 
ture (b) which has the greater NO distance and 
thereby stiffen the NN bond considerably. 

It was, therefore, necessary to calculate normal 
coordinates as a function of the variable fiz. The 
results of these calculations are summarized in 
Table IV, the indices being selected so that: 


vy; =2224 cm—, 0,=(L 4) unit (L") 1272. 
vo=1285cm-, Qe=(L™)eirit(L) oore. 


For the normal coordinates corresponding to 
each possible choice of fiz there are four possible 
combinations of signs for the du/0Q’s that must 
be considered. The values are: 


du /Qi= +244, du/dQ2=+111. 


The choice is restricted, however, by the fact 
that stretching the NO bond will favor structure 
(b) above, and, ther€fore, give the molecule a 
dipole moment directed to the left, whereas 
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Fic. 4. Hypothetical curves for dipole moment, y, as a 
function of internuclear distance, r. 


stretching the NN bond will favor structure (a) 
and give a dipole moment directed to the right. 
Substitution of the above values for du/0Q1 and 
0u/0Q2 with those from Table IV in Eq. (8) 
shows that this restriction is satisfied only if 
0u/0Q1 and du/dQ2 are both positive. Possible 
values of 0u/dr1 (NN bond), and du/dre (NO 
bond), are given in Table V as a function of fie. 
A positive contribution to the dipole moment on 
the part of the NN bond is defined as one which 
is directed from the central nitrogen to the end 
nitrogen. Similarly a positive contribution by 
the NO bond will be directed from nitrogen to 
oxygen. 

Table V also includes values for the dipole 
moments themselves, determined from the in- 
tensity of the bending mode. These signs are 
also uniquely determined. Since oxygen is more 
electronegative than nitrogen, the moment of the 
NO bond must be negative. It will be assumed 
that the molecular dipole moment is directed to- 
ward the nitrogen end and that it has the value” 
0.17<10—'8. These two conditions and the ob- 
served intensity determine the NO and NN bond 
moments using the hypothesis of additivity, 
though it should be remembered that in a reso- 
nant molecule any separation into individual 
bonds is somewhat artificial. 


Behavior of Bond Dipole Moments 


As might be expected, values of u/s and 
du/dr vary from about 0.1 to 2.0 times the charge 


a7 
2 J. H. Van Vleck, Electric and Magnetic Susceptibilities 
(Clarendon Press, Oxford, 1932), Chap. III. 
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of an electron. The slope, du/dr is in most cases 
greater than u/s, but the values for CH are 
consistent with a simple curve such as that 
shown solid in Fig. 4. The values for NN and 
perhaps NO seem to require a steeper form like 
one or the other of the dotted curves. The value 
of » for CH obtained from the 950 cm band is 
almost twice that obtained from the 1444 cm- 
or 995 cm bands. Apparently the electric 
moment behaves differently for this out-of-plane 
bending than for the in-plane modes because of 
some interaction with neighboring atoms. 


COMPARISON WITH OTHER DATA 


There are very few other measurements of 
infra-red intensities, and they are not directly 
comparable with those reported here. Bourgin™ 
obtained 0.86X10-" for du/dr in HCl, and 
Matheson" obtained 3.92 10- for CO. Foley" 
has given the ‘‘dipole moment” of the bending 
mode of HCN as 0.147 10-8, from which we 
calculate the CH bond moment as 0.57 10-3, 
u/s=0.54X10-. This is intermediate between 
our values for in-plane and out-of-plane modes 
in CoH,. 

The conventional value for the CH bond 
moment obtained from molecular moments" on 


®s £ 8 & oe 
i 
Aa 
ae 


$ 


x POINTS FROM 
EXPERIMENTAL CURVES 
OF ROLLEFSON ¢ HAVENS 


CONTRIBUTION TO (n#-1) IN UNITS OF 10°° 
4 : 
i=) 











; = 2 eS OU 
WAVE LENGTH IN MICRONS 
Fic. 5. Dispersion curve for ethylene calculated from 
intensity data and compared with experimental points 
from Rollefson and Havens. Dotted curves show limits of 
20 percent error. 


' iP. G. Bourgin, Phys. Rev. 29, 794 (1927); 32, 237 
1928). 

14. A. Matheson, Phys. Rev. 40, 813 (1932). 

15H. M. Foley, Phys. Rev. 69, 628 (1946). 

%C. P. Smyth, Dielectric Constant and Molecular 
Structure (Chemical Catalogue Company, New York, 





~ 1931), p. 117. 
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the basis of additivity is 0.4X10-, agreeing well 
with our in-plane moment. 

Measurements of infra-red dispersion also lead 
to values for bond dipole moments. In the case 


of methane Rollefson and Havens® concluded | 


that u=0.307 X 108, and du/dr = 40.562 X 10-™, 
in fair agreement with our in-plane moment and 
our values for du/dr. Rollefson and Havens have 
also very kindly furnished us with some un- 
published data on the index of refraction of 
ethylene, and we find that our values for the 
intensities reproduce their observed dispersion 
curve quite closely if the combination band at 
1890 cm is taken into account, as shown in 
Fig. 5. The agreement is well within an assumed 
maximum experimental error of 20 percent shown 
by the dotted curves except at the highest fre- 


TABLE VI. Atomic polarizations. 








Calculated from From visible index and 





Substance intensity data dielectric constant® 
C.H, 0.40 cc 0.39 cc 
N.O 0.55 cc 0.41 cc (0.46) 








®H. E. Watson and K. L. Ramaswamy, Proc. Roy. Soc. London, 
A156, 144 (1936). 
b E. Czerlinsky, Zeits. f. Physik 88, 515 (1934). 


quencies where contributions from electronic 
transitions should be taken into account. 

A final comparison can be made with the 
atomic polarization obtained from measurements 
of dielectric constant. Table VI shows the con- 
tribution of infra-red vibration bands to the 
molar polarization as calculated from our in- 
tensities and from the difference of extrapolated 
visible indices of refraction and measured di- 
electric constants.” Again the agreement is good. 
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On the Theory of Accidental Predissociation* 


I. KovAcs anp A. Bupé 
Department of Theoretical Physics, University of Budapest, Budapest, Hungary 


(Received November 11, 1946) 


Applying the results of quantum electrodynamics we investigate the case of perturbation 
between two discrete and one continuous term series. Even though no interaction between the 
first discrete and the continuous states exists, a considerable predissociation may arise for some 
spectrum lines in case of coincidence or close proximity between the levels of the two discrete 
terms belonging to the same rotational quantum number. The predissociation occurring in the 
first discrete term is caused by the second. This phenomenon is called accidental predissociation. 
Our results are in excellent agreement with the results obtained in the y-bands of NO molecule 
by Schmid and Gerd. 





HE phenomenon of accidental predissociation was first discovered by Coster, Brons, and van 
der Ziel’ in the so-called II positive band system of Ne. A short theoretical interpretation was 
subsequently given by Ittmann.? A similar effect has been observed recently by Schmid and Geré in 
the initial states of the y bands of NO.* The interpretation of these experimental data renders the 
extension and generalization of the theory necessary. 
The states of the system of the molecule and of the radiation field can be described according to 
Dirac‘ with a normalized wave function: 


(Ni, No, ---No, ++ +3 Gis Qa °° 3 15 fivie=t. (1) 
Here qi, g2, -:* are the space and spin coordinates of the electrons and nuclei, Ni, Ne, --- are the 
numbers of light quanta in the first, second, --- monochromatic plane wave, from which field, the 


radiation, is formed as a linear combination. Let us denote the wave functions of the molecular states 
without radiation field by ¢g:, where / is the totality of the quantum numbers related to the energy 
values. Since the ¢:’s (depending only on the space and spin coordinates of the particles) form an 
orthogonal, normalized complete system, we may write® 


¥(Ni, N», --N,, +S Gaels eee ;)=Da(M, No, ---N,, --+3f)oi(q1, q2, ee Z., |az|?= (2) 
Nol 


where the a: coefficients depend, apart from the time and the NV, quantum numbers of the radiation 
field also on the totality of the energy quantum numbers of the particles. 

Let us denote with H,, the operator corresponding to the separable part of the wave equation of 
the molecule and with G,, the energy operator of the radiation field, furthermore let H,,“ be the 
perturbation between the different levels within the molecule, which was neglected in setting up the 
function ¢, and let G,,“ be the small perturbing operator describing the interaction of the molecule 
and the radiation field. Then we have the time dependent wave equation as follows: 


h dy , 
[Hop + Gop + Hop +Gop YY = —— —. (3) 
2ni Ot 


* Lecture delivered by one of the authors (Kovacs) on September 25, 1946 in the Physics Section of the Hungarian 
Sciences Association. 

1C. Coster, F. Brons, and A. v. d. Ziel, Zeits. f. Physik 84, 304 (1933). 

2G. P. Ittmann, Naturwiss. 22, 118 (1934). 

8 R. Schmid and L. Geré, Math. u. Naturwiss. Anz. d. Ung. Akad. der Wiss. 62, 408 (1943). 

4P. A. M. Dirac, Proc. Roy. Soc. London A114, 243 (1927). 

5H. A. Kramers, Theorien des Aufbaues der Materie II. Quantentheorie des Elektrons und der Strahlung (1938), p. 462. 
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Substituting the expression for y from Eq. (2) into Eq. (3) and multiplying from the left by ¢* and 
integrating, we obtain the following expressions: 


fetta D1 aigidt =a,Ex, (4a) 
fercno pa aigidt=h(y,Nit+ voNeo+ ete +v,N,+ wifes *)ay, 2 (4b) 
fect. D1 aigidr =Yii Aridi, (4c) 


fotGu Li aigidt = Dip Li Ger [(N,+1)'ai(Mi, ---, Np+1, - +) 


+(N,)'ai(Ni, ree, N,-1, o° -)], (4d) 


h 
for (-—- — 1 aver )ar= ——dy, (4e) 
2771 Ot 271 


where E; are the unperturbed molecular energies belonging to the totality of the quantum numbers, 
>, hv,N, is the eigenvalue of the radiation field, Hi:= { ¢.*H., gidr is the perturbing matrix 
element of the non-radiating molecule and finally 


Gi” = f eG G2, -**)gidr. (4f) 


Here the operator G,, designates the part of G,,“ which depends on the space coordinates related 
to the vector potential of the pth plane wave.® With the help of the expressions mentioned above we 
have, from the wave equation, 


h 
—-—, = [E.th oe vpN, lar +>.1 Aha 


2a1 
+, D1 Gul (N +1)4ai(Mi, ---N,+1, -+-)+(N,)4a(Mi, ---N,—1,---)]- (5) 


For accidental predissociation, according to the experiments, we take into consideration the inter- 
action H4z“ between the two discrete term series denoted by A and B, furthermore the interaction 
Hzgc™ between the discrete states B and the continuous states C. We neglect the immediate inter- 
action between A and C. This means that the A state, in spite of the presence of the continuous state 
C, does not predissociate into C. The reason for this may be, e.g., that one of Kronig’s conditions on 
which the existence of the interaction depends is not satisfied. Further we make the supposition that 
an emission transition from the discrete term series A to a lower discrete term series denoted by D may 
exist. (See Fig. 1.) These assumptions mean that in Eq. (2) a; differs from zero only in the following 
cases: 


k=A, N,=N.=--:-=N,=-:-=0, 
k=B, N\=N2.=---=N,=--:=0, 
(6) 
k=C, NM,=N.=-::=N,=-:--=0, 
k=D, Ni =Ne2=:-+ =N,-1= Noy =:** =0, N,=1, 





*V. Weisskopf and E. Wigner, Zeits. f. Physik 63, 60 (1930). 
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~ A 3B Cc 
+e —EE 
Pi , 
Sy TaBLe I. 
SP os “ad 
ad ‘5 
Pe a 
J pert. pred. + ff oe Ea—Ep 
J f a an YY re fz lap, | 2dv fz ad,” | 2py 
| . a 0 cm 0.06 1.5 10-7 
/ —_o 2 cm™ 0.08 1.251077 
10 cm=! 0.24 3.2 10-8 
; 20 cm™ 0.5 10-8 
/ a 0 cm= 0.22 3.5 X10-8 
; Vee 5 cm= 0.34 1.9 10-8 
| 10 cm- 0.56 0.8 10-8 
' 20 cm>? 0.8 0.25 10-8 
J +i 
J-f L D 
3 
Fic. 1. 
furthermore 


Ea —Ep—he <hy, <Ex —Ept+he. 


With respect to the continuity of the state C we obtain from Eqs. (2) and (6) for y: 
& 3 p 1 2 p 
y=aa,(0, 0, -+-Q, -++3t)gatazn(0, 0, ---0, -++32) 


1 2 p 
+ fac(0, 0, ot . Es t)ec(EMdE+E, av (0, 0, --+,1,---3;Aegp. (7) 


Our problem is solved, if we succeed in computing explicitly the a; coefficients, because |a¢|* 
represents the probability of the realization of each state and if k=D, it gives the rate of decrease of 
the intensity itself. Substituting the expressions for y given by Eq. (7) into Eqs. (3), (4a), and (4e), 
respectively, Eq. (5) gives for the calculation of the a,’s the following system of linear differential 


equations: 
h 
“~—— =F4a4+Hapdept+ >, Gavan” (0---1---), (8a) 
1 
h 
~~ = Hpada +Exan+ { Hno(B)ac(E ME, (8b) 
1 
h 
——to(E) =Hcs(E)as+Eac(£), (8c) 
2771 
h ) 
——tip (0+ + +1- ++) =Gpaaa+[Ep thy, Jap (0- + +1+++). (8d) 


2a1 


We shall attempt to solve this system of equations mentioned above by the following assumptions: 
aa = Sae72tat : apg = Spe 2ret (9) 


where S4, Sz, and a are still to be determined. Substituting a4 into Eq. (8d), ag into Eq. (8c) with the 
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assumption that ap”) =a¢(E) =0, when ¢ is also 0 (that is, when only A and B states are excited) we 


obtain 
Ep 
1— exp( 2m] + r+ ialt) 








2ni 
ap” = SsGpa” -exp( -—“[Eothn}) ; 
Ep h 
i++ ia] 
h 
(10) 
2 
1-— exp(| 72 _ 2na +) 
(E) =Sp-Hoa(E) h ( -, 
ae =Op'fHics mm --— 1). 
hf 2x . h 
—|i=2- 2ra 
2miL h 


For the evaluation of S4, Sg, and the exponent a let us substitute Eqs. (9) and (10) into Eqs. (8a) 
and (8b). The summation over p in Eq. (8a), may be replaced, to a good approximation, by an 
integration over v between vp“’—e<v<vp4’+e.’ Considering that, according to the experimental 
results about accidental predissociation, the perturbed terms are hardly displaced from their original 
unperturbed places (in any case within the limits of experimental error), we may put |Gap |? before 
the integral sign with the value taken at y=vp4 (and not at v=yvp4’). Integrating in this way we 
have® 


h 
—i—Vape—****, ' (11a) 
ar 
where 
4r? 
Vav=—_ |Gan(vn*) |" (11b) 
1 


The physical significance V4p is the number of transitions from the A term into D per molecule per 
unit time (viz., the reciprocal of T4p, the average lifetime of A without predissociation). 
Proceeding in the same way, we obtain for the integration occuring in Eq. (8b)* 


h 
——Vgce-2™, (1 1c) 
4n 
where : 
An 
Veo=—~-|Hec(Es) |*. (11d) 
1 


Vec means the number of transitions from the B term into C without radiation per molecule per 
unit time (viz., the reciprocal of Tgc, the average lifetime of B without transition to a lower term). 
With these results we obtain from Eqs. (8a) and (8b) 


1 Vap 1 
s-E.+— a] +5, ‘-—H,ap=0, 
h 4n h 
(12) 
t 1 Vac 
Sa “Heat Sol Ext——a =(). 
1 


TT 


7 A’ is the perturbed value of the A term. 
® See reference 6, p. 63. 
*See H. A. Kramers, Theorien des Aufbaues der Materie. I. Die Grundlagen der Quantentheorie (1933), p. 219. 
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Equation (12) represents a system of homogeneous linear equations for the computation of S4, Sz. 
The condition of its solubility is that the determinant should vanish, v7z., 











1 Vap 1 

—-E4+——a —Hap 

h 4r h 
=(0 (13 
1 Vac 
—Hpa OO 
h h 4r 

Evaluating a, we have 
Vapt+Vae 2mrif Eat+Ep ” 
2ray y= ale) +—| Fignes?-a(0) | (14a) 


where 


A 1 47? €4? 2 H - 

oa mika aA) aia ae 
P 1 e4? 2 HW ? h 2 

a@=S|(—) +1 nen -(-«) 


(C2) sna) Taare) 0 


Vec—Vap 
ca? =Ep—E,; w=———. (14d) 


2,4 


+ | Hap} ‘= ew "| Has: — it (14b) 





Employing Eqs. (14a)—(14d), we obtain from Eq. (12) for the S’s (taking account of the conditions 
of normalization) . 





Sa=Se=(HI+(—e)'D; Ser=—Sae=(BL1- (1-0) 'T, (15a) 
where 
| Hap|? 
p= ee (15b) 
raw r+(—w) 


With the help of these equations we can determine the coefficients a, after some transformations in 
the exponent 
@ay=Saye*"*; Aan = Save", (16a) 


Qpi=Spie*"';  Ap2= Sp2e~**, (16b) 


10 More exactly we should have to write sign(e,?-w), but in our case, predissociation of the B state being complete, 
Vac> Vap, viz., w>0, we may neglect w. 
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2n1 
1 exp] —"E-a)t| 
2ni h 2ni 
aci(E) = Spi:—Hca(E) : -exp( -—"z1) > 
h 271 h 
—E-ay, 
h 
(16c) 
21 
1—exp| (“"z-ar)s] 
2ri h 2n1 
dc2(E) = Sg2:-—Hca(£) -exp( -- zt), 
h 2n1 h 
—E—ayz 
21 
1 ~exp(|="ep-+in,) -«:|t) 
2ni h 2n1 
ani” = Sa1-—Gpa” -exp( -—"[Ep-+H, 1) ; 
h 2ri h 
— rem 
(16d) 
2ri 
1 ~exp(|="ep-+in,) -«s[t) 
2Qri h 2n1 
apo” = S42-—Gpa” -exp( -—"[Ep+hn}), 
h 2r1 h 
—(Ep+hyv,) —az 
h 
where 
: 2rif Eat+Es ; 7 " 
ay = (S41? Vat Sar? Vac) +—| - —signes?- A(z) |, 
nL 
(17) 
aut Rat Be ? ; a 
a2 = 3(Sa2?Vap+Sz2" Vac) +—_ - > +signe,?- A(t) }. 
aL : 
In the expression (7) for ¥ the values of a, given in Eq. (16) are to be substituted. 
The |a,|*, that is, the probabilities of realization of each g; state, will be as follows: 
|@a1|?= Sai? exp(—[Sa1?Vant+Sai12Vec lt); |aa2|?=Sa2? exp(—[LS422Vap+Sp2?Vaec]t), (18a) 
|a1|?=Sp1° exp(—[Sar?Van+Sar°Vaclt); |ae2|?=Sp2* exp(—[Sa22Vant+Sp2*Vec}t). (18b) 
1 will give a detailed computation for one of the coefficients, 
; 2x rt” _, 1—2-cosyt-e~**+-e-* 
J laci(e) 4B = Sp," ~ f | Hen(E) |? dy 
Va y?+s? 
eae a dn Oe 
= Sp -—(1—e7*") = —exp(—LSa1’Vapt+Sai’ Vac Jt) |, 
; 2nr 2s Sar? eT ee, 
where 
2r E,t+Epg ; ’ ar 
E +signes?-a() |; $= 3(Say?>Vap+Sp1? Vac); dela 








172 I. KOVACS AND A. BUDO 


By a similar method we obtain for all the four coefficients 


Sar Vee 














Jlacice aE =~ 2Vapt+SeeV [1—exp(—[.S41?Vap+ Sai? Vac Jt) J; 
At” Vap+Sai1" Vee 
Spo? Vac aes 
Jaca) | "4B =— 2V7 1S 2y7 [1—exp(—[S42?Vap+ Spo? Vac It) |, 
Ao” Vap+SB2" V'Be 
(p) | 2g ~~ 1 Sar V. S 
p)| 2 = sian om 2 2/77 ° 
fz api | 2dv S°VaptSe'Veck exp(— [S412 Vap+ Sai? Vac |t) J; 
(18d) 
S42? Van 
fz oni iat 2 45 2) [1—exp(—[.S422?Vap+Sp2?V ec It) }. 
A2”Vap+Spo" Vee 
DISCUSSION OF RESULTS 
1. At the beginning of the time measurement (¢=0), we have: 
|@a1|?=Sa1?; |@a2|*=Sa2?, (19a) 
|@zi1|?=Sa1?; |a@p2|?=Sp2’, (19b) 
f lacs) |*4E=0; J lacs(E) |*4z=0, (19c) 
[Eo lan lde=0; [Xp laps |"dv=0, (194) 
and consequently 
¥i=SaigatSpigs; w2=Sa2gat Serge. (20) 


Disregarding the probability of the transition of the A term into D and of the B term into C 
(viz., Vap= Vac=0), Eqs. (20) means nothing else but the time-independent solution of the wave 
equation for the case of common perturbation between the A and B states. Namely in this case in 
Eq. (15b) the expression of p and A(z) will be: 


| Hap|? e4? 2 3 
pi ai)=((=) +1201") (21) 
| A(z) |? 2 





which are well-known results." 
2. In order to be able to compute the intensity decrease of the emitted light, we have to investigate 


our solutions for the case t= ©. Then we have: 














|aa1|?=0; |aa2|?=0, (21a) 
|az1|?=0; |az2|?=0, (21b) 
Spi? V. Spo” V, 
f tacs(&) |*4B8=——* Jaen) |WE=--— ate 
Sar?Vapt+Sa1?Vac Sa2? Vap+Sp2? Vac 
SarV. S42V. 
fx. | api | *dv= —— fz l|ap2™ | *dv= i ’ (21d) 
Sar?Vapt+Sai?Vec Sao? Vap+Sp22V ee 


1G, P. Ittmann, Zeits. f. Physik 71, 623 (1931). A. Budé and I. Kovacs, Zeits. f. Physik 109, 393 (1938). 
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Equations (21c) gives the percentage of the initial molecules which are transferred without radia- 
tion (predissociation) from the two perturbed states into C. As we may see, although at the beginning 
predissociation came only from one of the two unperturbed states, namely from B(Vzc0), now both 
of the two perturbed states are predissociated. 

Equations (21d) gives the percentage of the initial molecules which are transferred with radiation 
from the two perturbed states into the lower state D. Though at the beginning, transition was 
possible only from one of the two unperturbed states, namely from A(V4p#0), now each of the two 
perturbed states furnishes a spectrum line. 

Let us now consider what kind of numerical results can be deduced from our formulae for the 
decrease of the intensity. From other considerations it follows that under normal conditions we have 
Vap=108 (that is 74p=10-*), and in the case of a complete predissociation Vg¢=10'? (that is, 
Tsc=10-!). In the case of an accidental predissociation, according to the experimental results, we 
have to assume the H,4g matrix element to be so small that no displacements of the lines are ob- 
servable. Thus, in the first ease H4z should be 0.2 cm—", in the second case 0.1 cm! (both measured in 
wave numbers). Calculated for different term differences and using the above mentioned Hz values, 
we get the results given in Table I. 

From the data collected in the table it is obvious that the intensity of the lines corresponding to the 
perturbing state designated by the subscript 2 (corresponding to B) is negligible, while the decrease 
of the intensity of the lines arising from the perturbed state designated by the subscript 1 (corre- 
sponding to A) can take any values between 0.06-1.00 depending on the energy difference. 


COMPARISON WITH EXPERIMENT 


Schmid and Ger6é" photographed and analyzed the spectrum of the NO molecule between 1950- 
2700A with large dispersion and great intensity. They found that near the intersection of the v’=0 
level of the so-called 8-bands with that of the y-bands (using the (0,0), (0, 1), (0, 2), (0, 3), (0, 4) 
y-bands), there is no irregularity to be found, a behavior which can be explained by the small value 
of the matrix element of interaction. We find the same situation at the intersection of level v’=1 of 
the 6-bands with the level v’=0 of the y-bands at about J ~50. However the intensities of the lines 
belonging to K =48 and K =65, respectively, of the v’=1 and v’ =0 levels of the y-bands, are smaller 
than the intensities of the lines belonging to the neighboring quantum numbers; whereas while the 
wave numbers of these lines do not show any irregularities. These two lines correspond to levels just 
at the place where the v’ =2 and v’ =3 B-term series, which are already predissociated at this energy, 
cross the y-term series, leading to an accidental predissociation. Considering that even here the wave 
numbers of the lines exhibit no irregularity, we have to assume the matrix element to be very small, 
which seems to be a general property at these intersections corresponding to the case above mentioned 
(compare the crossing of v’=0 of y and v’ =0 of 8, and that of v’ =0 of y and v’ =1 of 8) This example 
is more convincing than that of Coster, since here the terms (the 6-terms) causing the accidental 
predissociation were in fact observed below the limits of their own predissociation (50,000-51,000 cm~'). 


® The work of Schmid and Geré quoted under reference 3 will be published later elsewhere. Furthermore see: K. Bal- 
lenegger, Mat. és Fiz. Lapok, Budapest 51 (1944). 
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I. Cation Arrangement in Spinels 


E. J. W. VERWEy AND E. L. HEILMANN* 
Natuurkundig Laboratorium der N. V. Philips’ Gloeilampenfabrieken, Eindhoven, Nederland 


(Received November 29, 1946) 


By means of a detailed x-ray study of a large number of oxides having a spinel structure, 
it has been established whether they crystallize in the ‘‘normal’’ spinel-type structure, XY20,, 
or in the “inversed”’ type of structure, Y(X Y)O,, discovered by Barth and Posnjak. Certain 
rules can be derived which are also applicable to more complicated spinels or to solid solutions 
of different spinel oxides. Madelung constants have been calculated for different spinel types in 
order to investigate their stability on the basis of an ionic concept of the chemical bond (Born’s 
lattice theory). The experimental results are only partially explainable on this basis. All 
exceptions to the electrostatic principles can be understood by assuming that Fe**, Ga**, In**, 
Zn?*, and Cd?* favor a configuration having the coordination number 4, because of a con- 
siderable covalent contribution to the chemical bond for that configuration. 





1. INTRODUCTION 


HERE are a great many oxides having the 
general formula XY2O, which crystallize 
with the same general crystal structure as the 
mineral spinel, MgAlI,O,y. Certain groups of 
spinels show remarkable electrical and magnetic 
properties. The spinels are, therefore, of con- 











Fic. 1. Arrangement and packing of atoms in unit cell of 
ideal spinel-type (u=0.375) containing 8[XY20,]. Large 
spheres represent oxygen, small black spheres represent 
fourfold coordination position and cross-hatched spheres 
sixfold coordination position. 


* Deported to a concentration camp in the eastern part of 
Germany (July 1944), from which he never returned. 


siderable interest from a technical point of view. 
They may be readily synthesized so that their 
properties, which are composition-sensitive, may 
be controlled within varying limits. 

The spinel structure is cubic with a large unit 
cell containing 8 X, 16 Y, and 32 oxygen atoms. 
The positions of the oxygen atoms are more or 
less fixed but the arrangement of the cations 
varies considerably within certain limits. The 





Fic. 2. Same as Fig. 1 but with unit cell shifted by one- 
half a cube edge. Small dashed-line cube shows NaCl-type 
structural motif with the cross-hatched ion in an octahedral 
interstice. 
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variations in distribution of cations in the X 
and Y positions cause marked changes in some 
of the physical properties. One of these is the 
electrical conductivity which will be considered 
in a following paper (II).. Equivalent lattice 
points may be occupied by different cations 
distributed statistically throughout the crystal. 
The permutations are often more complicated 
than those found in ordinary primary solid 
solutions. 

We have investigated a large number of these 
compounds during the past few years. The 
present paper describes a detailed study of the 
structure of the oxide spinels. It contains the 
crystallographic results for those compounds in 
which we have established the arrangement of 
the cations and a general discussion of the factors 
determining their arrangement. 


2. DISCUSSION OF THE SPINEL STRUCTURE 


In most oxidic structures the oxygen ions are 
appreciably larger than the metallic ions and 
the spinel structure can be approximated by a 
cubic close packing of O?- ions in which the 
Mg?+ and Al*+ ions occupy certain interstices. 
Each unit cell contains 8 [MgAlI.0,] and, there- 
fore, 32 O?- ions. This close packing contains 64 
interstices surrounded by 4 O?- ions (coordina- 
tion number 4, tetrahedral) and 32 interstices 
surrounded by 6 O?- ions (coordination number 6, 
octahedral). In the spinel structure, 8 of these 
tetrahedral holes and 16 of the octahedral holes 
are occupied. Unit cells of the “‘ideal’’ spinel 
structure are shown in Fig. 1 and Fig. 2; the 
only difference between them is that the cell 
unit has been shifted over one-half of the cube 
edge a in one of the three cubic cube edge direc- 
tions in the crystal. It is convenient to subdivide 
the unit cell in 8 octants. The structure then 
may be described as a face-centered compilation 
of A- and B-cubes, where the A-cubes comprise 
a group MgO, (Mg in 4-coordination) and the 
B-cubes a group Al,O, (Al in 6-coordination).! 

In all spinel-like substances the parameter u 
has a value in the neighborhood of 0.375. For 
u=0.375, the case corresponding to Figs. 1-3, 
the arrangement of the O?- ions equals exactly a 
cubic close packing. In actual spinel lattices this 


1E. J. W. Verwey, Zeits. f. Krist. 91A, 65 (1935). 















































Fic. 3. Two octants of a unit cell. Spinel structure is a 
three-dimensional packing of A- and B-cubes in a checker- 
board arrangement. 


ideal pattern is slightly deformed, usually, corre- 
sponding to u>0.375, in such a way that the 
oxygen tetrahedron in the A-cubes is somewhat 
expanded and the oxygen tetrahedron in the 
B-cubes is slightly contracted. Accordingly, in 
that case, the octahedrons formed by the 6 O?- 
ions directly surrounding the positive ions in 
these B-cubes deviate somewhat from regular 
octahedrons. 

In simple spinels the valencies of the cations 
may be either 


X=2+, Y=3+, 


X=4+, Y=2+. 


The third possibility, X=6+, Y=1+, is very 
rare. In the present work, separate examples of 
the two common possibilities and their combina- 
tions in more complicated spinels are studied. 
Originally, in accordance with classical prin- 
ciples of crystallography, it was assumed that in 
the spinels the X atoms or ions occupied the 
eightfold position (the tetrahedral holes of the 
A-cubes indicated in Fig. 3) and the Y atoms or 
ions were in the 16-fold positions (the octahedral 
interstices of the B-cubes). However, Barth and 
Posnjak? pointed out a second possibility of 
distributing the cations while retaining cubic 
symmetry; viz., one in which the tetrahedral 
interstices are occupied by half of the Y ions, 
the octahedral holes by the X ions, and the other 
half of the Y ions are distributed at random at 
these lattice points. For the sake of convenience, 
the first arrangement will be indicated as charac- 
teristic for “normal spinels,” the Barth and 
Posnjak arrangement as that of “‘inversed spinels.” 
A number of examples for both groups have 


2 T. F. W. Barth and E. Posnjak, Zeits. f. Krist. 82, 325 
(1932). 
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already been given by Barth and Posnjak, who 
based their conclusions upon a detailed dis- 
cussion of the intensities of x-ray diffraction 
pictures. In addition to these two arrangements, 
there is the possibility of some intermediate 
arrangement with an averaged distribution of all 
ions about all spinel cation positions. As a matter 
of fact, this contingency constitutes in effect an 
infinite series of possibilities, as this type of 
distribution will generally occur in such a way 
that certain ions show a definite preference for 
' one type of interstice according to a certain 
distribution factor. Machatschki* has described 
such an intermediate case (MgGa,0,). 


3. CONSIDERATION OF THE STABILITY OF DIF- 
FERENT ARRANGEMENTS BASED ON CALCU- 
LATIONS OF MADELUNG POTENTIALS 


Before passing to an extension of the experi- 
mental facts and an attempt to find general 
rules as to the most stable arrangement in a 
given case as a function of the properties of the 
participating ions, it is useful to consider first 
what might be expected from a purely heteropolar 
point of view. 

The chemical bonding in these spinel oxides 
can be treated to a first-order approximation as 
an ionic bond, in which the state of affairs is 
described by an equilibrium of Coulomb forces 
and Born repulsion forces. Comparing states 
with the same crystal structure, the former forces 
are described by the Madelung potential; the 
latter are roughly determined by the size of the 
different ions. The ionic radii will not vary over 
a very wide range, because the spinel structure 
is stable only if the cations are rather medium 
sized and, in addition, the radii of the different 
ionic species in the same compound do not 
differ too much. The ionic dimensions will have 
an influence upon the value of the parameter 1; 
large values of u can be expected if the tetra- 
hedral interstices are occupied by comparatively 
large ions. Values for this parameter, however, 
have been determined in only a few cases (where 
they were usually about 0.380). Accurate values 
are ‘rather difficult to obtain and only very few 
are available. 

Madelung summation constants have been 
calculated for spinels with the following values 


3 F, Machatschki, Zeits. f. Krist. 82, 348 (1932). 
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of the valencies, where p and g are the valencies 
of the ions at tetrahedral and octahedral inter- 
stices, respectively : 


(a) p=4, q=2, 
(b) p=2, g=3, 
(c) p=3, g=23. 


Case (a) will then give the Madelung constant 
of a normal 4-2 spinel, X**Y,.?*O,, and case (b) 
that of a normal 2-3 spinel, X?*Y,.*tO,. In the 
inversed 4-2 spinel Y?*(X**Y?*)QO,, the 16-fold 
position contains an averaged distribution of 
4-valent and 2-valent ions and the corresponding 
Madelung constant will approximately be given 
by case (b) as the valency of the ions in the 
16-fold position is on the average 3. In the 
inversed 2-3 spinel Y**(X?*+Y**)Q,, the valency 
between the brackets is on the average 2} and 
the Madelung constant will be approximated 
by case (c). 

The calculations have been carried out for 
different values of « in order to establish the 
effect of this parameter upon the stability. 

The calculations have been done by a com- 
bination of the Ewald‘ and the Evjen® methods, 
for the following reasons. As a consequence of 
the large number of ions per unit cell the calcu- 
lations are rather laborious, especially in those 
cases where the parameter u differs from 0.375 
and the oxygen ionic lattice shows small devia- 
tions from the ideal face-centered cubic com- 
pilation. 

Evjen’s method, a direct summation over a 
restricted number of unit cells where electro- 
neutrality is maintained by counting the ions in 
the outside planes 3, etc., appears to lead to 
results accurate within 1 percent when the 
summation is carried out over a cube of 8 unit 
cells, with the ion under consideration in the 
corner common to all 8 cells (this has been 
checked by extending the calculations in one 
case over 27 and 64 cells). Accordingly, this 
method has been used for a calculation of the 
Madelung constant M for u=0.375, 0.380, 0.385, 
and 0.390. In some cases, including the present 
one, however, the Evjen method leads to M 
values which contain a constant error due to 


4P, P. Ewald, Ann. d. Physik 64, 253 (1921); Zeits. f. 
Krist. 56, 129 (1921). 
5H. M. Evjen, Phys. Rev. 39, 675 (1932). 
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the circumstance that the outer part of the 
crystal fraction cut out of the crystal in this way 
gives rise to an electrical dipole layer. As this 
constant error will be independent of wu, it is 
sufficient to eliminate it by calculating for the 
ideal parameter value u=0.375, the Madelung 
constant according to an alternative method 
avoiding this constant error. For this calculation 
Ewald’s method has been chosen. According to 
this method the summation over the infinitely 
large lattice is carried out by transforming the 
infinite series into the elliptical theta-functions, 
so that again a rapid convergency is obtained. 
Calculating in this way the electrical potential 
at a tetrahedral lattice point, an octahedral 
lattice point, and an O-point, we find for the total 
Coulomb energy per ‘‘ molecule’ XY2O,: 


V=- 4 F (e/a) (DM tetra. +2qMoct. +8Mo) 

= —M- (e/a) 
where a=the length of the unit cell. The data 
obtained for the Madelung constant M are given 
in Table I. Hence, on the basis of the ionic picture 
of spinels, assuming that the cube edge of the 
unit cell is not altered by the inversion (we will 
see later on that this is not entirely true, but 
sufficiently true for our following conclusions), 
and also assuming that the Born repulsion 
potential is a constant fraction, amounting to 
about 10 percent of the Coulomb potential, we 
should arrive at the following results concerning 
the stability of different cation arrangements. 

1. The normal 2-3 spinel (case (b)) is more stable than 
the inversed 2-3 spinel (case (c)), but the inversed 4-2 
spinel (again case (b)) is more stable than the normal 4-2 
spinel (case (a)). 

2. This conclusion is not altered by the influence of the 


parameter, as the M values do not cross within the range of 
practically occurring u values. 


The first conclusion can also roughly be ex- 
plained on the basis of the principle of maximal 
charge neutralization ; in both cases we find that 
the most stable arrangement is that in which 
the cations with the highest charge occupy the 
interstices with the highest coordination number. 

Similarly it can easily be understood that the 
Madelung potential of a normal 4-2 spinel is 
increased by a decrease of u and that of a normal 
2-3 spinel by an increase of u, as both variations 
of u correspond to a contraction of the oxygen 
ions around the cations with the highest charges. 











TABLE I. 

(a) (b) (c) 
p=4; p=2; p=3; 
q=2 q=3 q=23 

u=0.375 142.1 150.3 143.6 
=0.380 139.9 153.1 143.8 
=0.385 137.7 155.7 144.1 
=0.390 135.6 158.2 144.3 








It is also understood that the intermediate case 
(p=3, g=234) is less sensitive to variations of the 
parameter. 

In the following section the extent of the 
agreement between the predictions given above 
concerning the relative stability of different 
arrangements and the cation distributions ac- 
tually found will be shown. 


4. X-RAY INVESTIGATION OF VARIOUS SPINELS 


A large number of spinel oxides were prepared 
by mixing the pure oxides (‘‘pro-analysis”’ 
quality) in the proper ratio, pressing the mixture 
into bars and heating the latter for several 
hours in an electric furnace at a constant tem- 
perature. The sintering temperature varied with 
the nature of the constituent oxides; for easily 
sintering samples it was mostly 1200°C, for 
oxides sintering with greater difficulty it was 
1300°C, or even 1400°C (as in the case of 
MgCr.0,). Only in the instance last mentioned 
was it necessary to carry out the heating process 
in an atmosphere of oxygen; in the other cases 
air was sufficient to prevent reduction. The 
samples with Fe*+ form an exception to this 
rule and require special precautions to obtain 
the iron completely in the desired valency state. 

The determination of the cation arrangement 
and, consequently, also of the detailed pattern 
of the oxygen ions (the latter described by the 
parameter “) was made in the usual way with 
the aid of x-ray powder diagrams.* The in- 
tensities of a great number of diffraction lines 
were measured and the observed values were 
compared with the calculated reflection intensi- 
ties** for different arrangements and different 

* We wish to thank Dr. J. F. H. Custers and Mr. C. G. J. 
Jansen, both from this Laboratory, for the preparation of 
these x-ray pictures under carefully controlled, checked 
conditions. 

** The intensities were measured visually by comparison 
with an intensity scale. In the calculations the usual angle- 


dependent factors were taken into account, including 
absorption. 
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TABLE II. Intensities. 











Plane Cu Cuo.sZno.s Cuo.«Zno.6 Zn 
(111) 5 4 3 1- 
(220) 15 20 20 30 
(311) 100 100 160 100 








values of u. In several cases the difference in the 
scattering power of the constituent metal ions 
for arbitrarily chosen wave-length of x-radiation 
was too small to allow a decision between differ- 
ent arrangements. In that case it was often 
possible to arrive at the desired result by a 
special choice of the irradiating wave-length. 
The scattering power of an element is compara- 
tively low for the characteristic Ka-radiation of 
the metal having an atomic number one higher 
than the scattering metal atom, as a consequence 
of the proximity of the absorption edge. Thus, 
the scattering power of Fe is reduced by about 
20 percent for the same values of sin@/A when Co 
Ka-radiation is used instead of, for instance, 
Cu Ka-radiation. Hence, we must replace the 
usual theoretical scattering factor by the experi- 
mental values, which in this case have been 
determined by Bradley and Hope.® In this way it 
was possible, for instance, to determine the 
arrangement in samples containing both Cu and 
Fe which normally would have been impossible 
to determine. 

The following substances were investigated: 
MgCr2Ou, NiCr2Q,, CdCr20,, CuFe.0,, ZnFe20,, 
Cu,Zn;_zFe20u, CdFe,0,, TiZn2O,, SnZn2O,, 
SnCo.0,, Fe?+Fe*+AlO,, LiAl;Os, and a large 
number of solid solutions containing these and 
other spinels. 

The following are the results obtained with the 
simpler spinels. 

MgCr.0,4. Radiation Cu Ka. Lattice dimen- 
sion a=8.312A. Intensities exclude the arrange- 
ment Crs(MgsCrs)Oz2 and point definitely to 
MgsCrieQOz2, hence all Mg in eightfold position 
(tetrahedral holes), and all Cr in 16-fold position 
(octahedral holes). Comparison of experimental 
and calculated intensities does not permit very 
accurate conclusions as to the value of u; the 
calculations were made with differences of 0.5/60. 
The best agreement is found for u = 23/60 though 


6 A. J. Bradley and R. A. H. Hope, Proc. Roy. Soc. 136, 
283 (1932). 
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also u=23.5/60 cannot be excluded. Hence 
u=0.38;+0.01. 

NiCr,O,. Radiation Cr Ka. a=8.299A. For the 
scattering power of Cr for Cr Ka a correction was 
made in view of the vicinity of the absorption 
edge. Arrangement definitely NigCri¢O3e. Intensi- 
ties do not permit a decision between u = 23/60 
and u=23.5/60, hence u=0.38, to 0.39+0.01. 

CdCr,O0,. Radiation Cu Ka. a=8.567A. Ar- 
rangement CdsgCrigQ32. Best u value 23/60 
though 23.5/60 cannot be excluded. Hence 
u=0.38;+0.01. 

CuFe,O.. Radiation Co Ka. a=8.37A. Very 
good agreement for inversed spinel arrangement : 
Fes(CugFes)O32 with u = 0.38 )+0.005. 

ZnFe,0,4. Radiation Co Ka. a=8.416A. Very 
good agreement for normal spinel positions: 
ZngFeyg¢0 32 with u=0.385+0.005. 

Cuo.sZno.sFe2O4. Radiation Co Ka. a=8.388A. 
Arrangement Zn4Fe4(CusFe12)O32 with u =0.380 
+0.005. 

Cuo.sZno.6Fe2O.. a@=8.402A. Arrangement 
Zn4.sFe3.2(Cug.2Fe12.3)O32 with «=0.380+0.005. 

It is not possible to distinguish between Cu 
and Zn. However, the gradual transition between 
the inversed spinel (copper) and the normal 
spinel (zinc), with each Zn and Cu keeping its 
own position, can very well be followed in the 
gradual change of the intensities as one goes 
from CuFe,0, to ZnFe,O,. For example, the 
relative increase of the intensity of (220) or 
decrease of (111), in comparison with (311), 
which was a strong reflection arbitrarily set equal 
to 100 in all films, is shown in Table II. 

CdFe,O,. Arrangement CdsgFe1¢Q0 32. 

TiZn,O4. Radiation Cu Ka. a=8.445A. Ar- 
rangement Zng(TisZng)Oz2 with u=0.380+0.005. 

SnZn.O,. Arrangement Zng(SngZng)Ozo. 

SnCo,0,. Arrangement Cog(SngCog)O32. 

Fe*+Fe*+AlO,. Radiation Co Ka. 

Best agreement for the arrangement 
Fes?*(AlsFeg*+)O39 or Fes*+(Als*+Fes?*)O30, with 
u=0.385+0.01, though small deviations remain 
which can be explained by a small percentage of 
Al in eightfold position (less than 10 percent). 
The presence of a weak (111) reflection might 
even suggest that Fe?+ is at the 16-fold position, 
as for small reflection angles the scattering 
power of Fe?+ will be appreciably larger than 
that of Fe**. 
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LiAl,;Og. The x-ray data used in the investiga- 
tion of this substance were those found in the 
literature.’ Of the two possible arrangements, 
Alg(LigAl12)O32 or (LigAl4) Al 16030, the intensities 
are definitely in favor of the former, though the 
best agreement.is obtained by assuming that the 
actual distribution is an intermediate one, with 
a small contribution of the second, 4=0.385 
+0.01. Hence, the most probable arrangement 
approaches the completely random distribu- 
tion which would correspond to (Lii 1/3Al¢ 2s) 
(Lie 2/3Ali1s 1/3) Op. As this substance is the only 
one investigated of the interesting type in which 
Lit +Al** together take the role of two divalent 
ions, we merely give the results for this com- 
pound; in the following sections, which are 
concerned only with 2-3 and 4-2 spinels, it will 
be omitted from further consideration. 


5. CONSIDERATION OF LATTICE DIMENSION; 
DISCUSSION 


Including the results of Barth and Posnjak* 
(who found normal spinels with a number of Al 
spinels, inversed spinels for Mg2TiOu, Fe2TiOz, 
MgFe2Ox,,and MgIn.O,) and that of Machatschki 
(MgGa,O, intermediate case), we arrive at the 
following conclusions: 


(a) 4-2 spinels are always found in the inversed arrange- 
ment, in complete accordance with the lower Coulomb 
potential thereby realized. 

(b) 2-3 spinels are often found in the normal arrange- 
ment, again in accordance with expectation on the basis of 
electrostatic considerations. Thus, all Al and Cr spinels 
investigated obey this electrostatic rule. The same holds 
for ZnFezO,4 and CdFe:O,. Obvious exceptions, however, 
are a number of other Fe spinels (MgFe2O., CuFe.0,) and 
also MgInzO4, MgGa2O, (partially inversed). 


In connection with these conclusions, it would 
be interesting to have more information as to 
the behavior of other Fe spinels, such as those 
containing Mn, Fe, Co, Ni as the divalent ion. 
Here, however, the difference in scattering power 
isso small that no information can be obtained 
from the x-ray intensities. We have investigated, 
therefore, whether there are regularities in the 
values of a, the cube edge of the unit cell, for the 
different spinels. This is shown in Table III, 
where most data concerning @ were derived from 
the literature; in those cases for which more 


7E. Kordes, Zeits. f. Krist. 91, 193 (1935). 
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than one value was found we used an average 
value. 

Comparing all Al spinels, all Cr spinels, and 
those Fe spinels for which a normal arrangement 
was established in the present paper (Zn and Cd), 
it is seen that the lattice dimensions of these 
spinels run completely parallel to their composi- 
tion and, apart from small deviations, go in the 
sequence to be expected from the divalent ionic 
radii. One finds an increase of 0.235A (average 
value maximal deviation only 0.015A) in a as 
one goes from Al to Cr spinels, and an increase 
of 0.12A as one goes from Cr to normal Fe spinels. 
This increase, practically independent of the 
nature of the divalent ion, can easily be under- 
stood on the basis of the different sizes of the 
Al**, Cr*+, and Fe** ions, and it is in agreement 
with the gradual increase of the ionic radii 
established from more simple oxide structures. 

Making the same comparison, but now with 
the other Fe spinels, we find the same parallelism 
and again a constant increase of a between Cr 
and inversed Fe spinels. Now, however, this 
increase is only 0.04 to 0.06A and, accordingly, 
on the average is 0.07A less than that found for 
the two normal Fe spinels (Zn and Cd). Since 
the Fe spinels of Cu and Mg were found to be 
inversed spinels, this may tentatively be con- 
sidered as an indication that the Fe spinels of 
Mn, Fe?*, Co, and Ni are also inversed spinels; 
simultaneously one can conclude that the Al and 
Cr spinels for which the normal arrangement was 
not established by x-ray intensity considerations 
are actually all normal spinels. 


TABLE III. Cube edge a of unit cell for a number of 2-3 
spinels X Y 20x. 








yst = 


ryit = 


AB+ 
0.57 


Crt 
0.64 





a 


X2+ rx2+ a a 





Ni2t 
Cu?t 
Zn2+ 
Mg? 
Co?+ 
Fe?+ 
Mn?2t+ 
Cd2+ 


8.36 
8.37 
8.42 
8.36 
8.36 
8.39 
(8.55)** 
8.69 


8.30 


8.30 
8.31 
8.32 
8.34 
8.49 
8.57 


0.78 
* 


0.83 
0.78 
0.82 
0.83 
0.91 
0.97 


0.05 
(0.06) 
0.12 
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* The value of Tout =1.01 found in many tables in the literature is 
obviously incorrect. : : 

** The values for MnFe20, in the literature vary considerably, 
obviously as a consequence of differences in the valency states of Mn 
and perhaps also Fe. 
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TABLE IV. 








Solid solutions of: Formula 





Fe3+(Fe?+Fe%t)O,4 and all a Fe3+(Fe2+Al)O4 

Fe (Cu Fe )Osand Zn (FeFe)Ou (Fe1/2Zn1/2)(Fes;2Curs2)O% 
Fe (Mg Fe )Osand Mg (CrCr)O4 Fe (Mg Cr)O« 

Zn (Fe Fe )Osand Zn (CrCr)O« Zn (Fe Cr)O« 

Zn (Ti Zn )Osand Fe (MgFe)O« Zn (FeMgi2Ti12)Os 
Mg (Ti Mg )Osand Fe (MgFe)O. Fe (Tiij2Mgaj2)O4 


Zn (Fe Fe )Osand{Re Fierys)0s (Znu2Fer2)Fei 960s 








Hence, we arrive at the conclusion that, 
contrary to electrostatic rules, all Fe spinels are 
inversed spinels, except those of Zn and Cd. 
The comparison of lattice dimensions is obviously 
not so conclusive as a direct derivation from the 
intensities, but our conclusions will be seen to 
be in complete accordance with other physical 
properties of the Fe spinels. 

The shrinkage of about 0.07A found above for 
an inversed spinel in comparison with the corre- 
sponding normal arrangement must obviously 
be explained by a more dense structure in the 
former case. Actually, the arrangement which is 
electrostatically most stable for a 2-3 spinel, 
where the comparatively large divalent ions go 
into the smallest interstices, is not the most 
favorable one from a geometrical point of view, 
and a certain contraction of the whole lattice 
can reasonably be expected if, for some reason, 
the inversed spinel becomes the most stable con- 
figuration (the corresponding contraction is itself 
a factor in favor of the inversed configuration, 
but it is sc small that its effect is negligible in 
comparison with the difference between the two 
Madelung potentials). 

Since it is thus made plausible that all Fe 
spinels, except Zn and Cd, have the inversed 
cation spinel arrangement, the explanation of 
this abnormal behavior must obviously be found 
in the circumstance that these structures cannot 
be completely understood from a purely ionic 
point of view. Comparatively important homo- 
polar contributions to the bond must be con- 
sidered. Apparently the Fe*+ ion shows a certain 
preference for such a homopolar contribution if 
surrounded by 4 O? ions, and the same must be 
true for In*+ and Ga**. The fact that ZnFe:O, 
and CdFe2O, show again a normal behavior can 
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then easily be explained, as the Fe** will find in 
Cd?* and Zn** a competitor for the tetrahedral 
interstices. The preference of Zn*+ and Cd?* for 
a tetrahedral surrounding is well known; e.g., 
in ZnO, ZnS, Zn2SiOu, CdS, CdSe; and has been 
attributed to homopolar contributions to the 
chemical bond.’ Zn?+ and Cd?* have a filled d- 
band, and the distance between the d-band and 
the unoccupied s-band is only small; hence there 
is a possibility of hybridization of 3 d-electrons 
and 1 s-electron, forming 4 bonds in tetrahedral 
directions (sp* and d*s form a hybrid with 4 
equivalent tetrahedral bonds’). The samé expla- 
nation would probably hold for In*+ and Ga*+ 
which have the same electronic configuration as 
Cd?+ and Zn?*. 

The general principles underlying the building 
up of oxidic spinels can thus be summarized as 
follows: 


1. Trivalent and quadrivalent ions, in accordance with 
electrostatic rules, preferably occupy the octahedral inter- 
stices of the 16-fold position. 

2. Exceptions are, however, certain trivalent ions: Fe**, 
In*+, Ga’+ which show a preference for the tetrahedral 
interstices of the eightfold position. 

3. Some divalent ions, Zn?+ and Cd?*, show a special 
preference for occupying the tetrahedral position, which 
can counteract the tendency mentioned under 2. 


These general principles appear to hold their 
validity in more complicated spinels. The appli- 
cation of these rules may be illustrated with the 
aid of a few examples given in Table IV.* 

The atomic configuration in all these cases has 
been confirmed by an analysis of the intensities 
of the x-ray diagram (for the first two cases this 
has already been given in Section 4 of the present 
paper) and other means of investigation. 

The crystallographic points of view developed 
here will be used in a series of subsequent papers 
in which a study of various physical properties 
of many complicated spinels will be reported. 


8 F. A. Kréger, Thesis, Amsterdam, p. 104 (1940). 

®G, E. Kimball, J. Chem. Phys. 8, 188 (1940). 

* In this and the following paper we use a notation which 
is slightly different from that in mineralogical work. 
Fe(MgFe)O,, for instance, expresses that the octahedral 
position contains together 1 Mg and 1 Fe for each Fe in the 
tetrahedral position. 
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The relations between the electronic conductivity of certain spinels and the arrangement of 
the cations in the crystal structure (see preceding paper) are studied. Several arguments favor 
the assumption that Fe;O, contains both divalent and trivalent iron in the 16-fold position. 
The transition point in the neighborhood of liquid-air temperature is probably associated 
with an increased degree of order at low temperature in the distribution of the 8 electrons 
between the 16 Fe-lattice points per unit cell. The considerably increased conductivity below 
the transition points shows tetragonal anisotropy when the crystal is cooled in a magnetic 
field. The possible distributions of the electrons in the crystal at low temperature are dis- 
cussed. In more complicated spinels, containing other metal atoms as well as iron in both the 
divalent and the trivalent state, the electronic interchange is more or less inhibited by the 
foreign metal atoms. The higher values of their resistance in comparison to that of FesO,4 can 
be roughly described by an increased activation energy. The investigation of a number of 
substances with different arrangements of the cations shows that the activation energy (and 
therefore the electrical resistance) is lowest for those cases in which the electrons can travel, as in 
Fe;0,4, along the Fe of the 16-fold position. 





1. CATION ARRANGEMENT OF Fe;0, sudden decrease of the conductivity was ob- 


served which was explained in terms of an order- 
disorder phenomenon of the Fe** and Fe** in- 
volved in the electronic conduction.* 


HE simplest compound for which inter- 

esting relations appear to occur between 
electronic conductivity and the arrangement of 
the cations is Fe;Ox. 

In paper I of this series,! we classified Fe;O,4 as 
an inversed Fe spinel Fe*(Fe*+Fe*+)QOx,, to- 
gether with the Fe spinels of Mn?*, Co, Ni, Cu, 
and Mg, as contrasted with the ‘normal’ spinels 
ZnFe2O, and CdFe2QO,. This classification could 
be made plausible on the basis of a detailed con- 
sideration of the values of dimensions of the unit 
cell of a large number of spinels. Several years 
ago, however, this “‘inversed’”’ arrangement had 
been already proposed by Verwey and de Boer? on 
the basis of the electrical properties of FegO4. The 
high electronic conductivity, ¢ = 200 ohm cm-, 
in comparison with that of similar compounds 
(Mn;3QOx,, y-Fe2O3) which have a very low con- 
ductivity, was explained on the basis of the 
continuous interchange of the electrons between 
Fe*+ and Fe** at the 16-fold position. This picture 


a. The Value of the Parameter u 


As an argument in favor of this inverse 
arrangement? the value of the parameter uw had 
already been used in connection with the ionic 
radii of the Fe?* and Fe** ions. This argument can 
now be further substantiated by the Madelung- 


potentials calculated in reference 1. This 
Madelung constant was given there as a function 
of u, and it is also possible to calculate the Born 
repulsion between the different iron ions and the 
surrounding O?- ions as a function of this 
parameter. The total lattice potential (Coulomb 
+ Born potential) can thus be calculated, for a 
given (i.e., experimental) value of the lattice 
constant a, and the value of uw for which this 
total potential is a minimum can be determined. 

For the repulsive potential we used the Born 


was corroborated by the behavior of Fe;O, at low 
temperatures, where below a transition point a 


1E. J. W. Verwey and E. L. Heilmann, J. Chem. Phys. 
15, 174 (1947). 

*E. J. W. Verwey and J. H. de Boer, Rec. trav. chim. 
55, 531 (1936). 


Mayer function 
Ve=bexp {(r+r_—r)/p} 
with for Fe?+ r, =0.83, for Fe** r, =0.67, and fora 


( sd 7h W. Verwey and P. W. Haayman, Physica 9, 979 
1941). 
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Fic. 1. Unit cell showing possible arrangement at low 
temperature of Fe** (large spheres) and Fe** (small spheres) 
in 16-fold octahedrally-coordinated position. All 32 oxygens 
and eight Fe** in eightfold tetrahedrally-coordinated posi- 
tion are omitted for sake of clarity. Projection of the atoms 
in the 16-fold position onto the back zy-plane and bottom 
xy-plane are shown at the left. The numbers indicate the 
fraction of a above these planes, respectively. The small 
black dots represent Fe**+ and open circles Fe**. 


mixture of Fe*+ and Fe?* the average value 
r,=0.75 ; for the oxygen ions r_ = 1.40; p=0.345A, 
the usual constant. Substituting u=%?+46, we 
obtain the total repulsive potential per molecule 
XY.0., on the assumption that the oxygen- 
oxygen repulsion forms only a small contribution 
to the potential and can be neglected as a first 
approximation : 


Ve =4b-exp[ {2.07 —av3($+4)} /0.345 ] 
+12b-exp[ {2.15 —a(4—6)} /0.345 ]. 


This function has a minimum value for 
5=0.002 or u=0.377. The attractive potential, 


however, 
Ve=—Me'/a 


is also a function of u or 6, as can be seen from 
Table I in reference 1. After eliminating the 
constant 6 with the aid of the equilibrium con- 
dition (dV /da) =0 for the experimental value of 
a, the total potential V=Vc+Ver appears to 
have its minimum for u=0.379. The experi- 
mental value found by Claassen‘ is u=0.379, 
from the x-ray intensities determined with the 
aid of an ionization chamber for a large crystal of 
Fe;0,4. For the normal spinel Fe?*(Fe**+Fe?+)O, 


4 A, Claassen, Proc. Phys. Soc. 38, 482 (1925-26). 


we would have calculated a much larger value of 
u, in the neighborhood of ~=0.390 (here the 
mutual repulsion of the oxygens cannot be neg- 
lected, owing to the much larger deformation of 
the close packing of the oxygen ions correspond- 
ing to the larger value of’ u). The calculation 
shows that the larger values of u to be expected 
for normal 2-3 spinels (and actually found for all 
Al spinels, etc.) are mainly caused by the 
Coulomb potential, the absolute value of which 
increases considerably with increasing value of w. 


b. Further Experiments on the Conductivity 


Assuming for FesO, the inversed spinel arrange- 
ment, the distribution of 8 Fe?*+ and 8 Fe** about 
the 16 lattice points of the octahedral position 
can be described as a distribution of 8 electrons 
about 16 Fe** ions. Above the transition tempera- 
ture this will approach a statistical distribution, 
because of the continual electronic interchange. 
Below this temperature, however, these con- 
ducting electrons will tend to arrange themselves 
according to some regular pattern. Without 
resorting to a larger unit cell this pattern of order 
can no longer be of cubic symmetry. The most 
probable pattern is one having a tetragonal ceil 
and is shown in Fig. 1. 

Hence the pattern of order can be described as 
being built up by continuous rows, consisting 
alternately of either Fe** or Fe** only, in the 
xy-plane. The Fe?* rows are all parallel to each 
other, as are the Fe**, but in a direction perpen- 
dicular to the Fe** rows. If the Fe?* ions in the 
plane z=% occupy the lattice points 4, 4; %, 3; 
5,2; 4, 3; the Fe? ions in the plane z=% occupy 
the points §, 3; 3) $3 8 83 8 & etc. 

If this description of the order-disorder transi- 
tion at low temperature is correct, one should 
expect that the electronic conduction below the 
transition temperature would be different in 
different directions. This will be measurable only 
if the direction of the tetragonal axis can be fixed. 
Ching Hsien Li’ has found that a single crystal of. 
Fe;0, becomes magnetically anisotropic at low 
temperature. Hence it seems possible to fix the 
tetragonal axis with the aid of a magnetic field. 
We did not have a single crystal of Fe;O, of 
sufficient purity to show the transition at low 


’ Ching Hsien Li, Phys. Rev. 40, 1012 (1932). 
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temperatures. In a previous paper*® we have 
shown that the transition phenomenon is entirely 
suppressed by a small amount of oxygen in 
excess of stoichiometric proportions; small 
amounts of impurities dissolved in the lattice 
have the same detrimental effect as a deviation 
from the stoichiometric composition. We had to 
use, therefore, sintered bars of pure Fe;O0u,, 
consisting of agglomerates of small crystals 
oriented in all directions, where the expected 
effect will be reduced considerably. The following 
experiment, therefore, has only qualitative 
significance. 

A sintered bar was cooled in a magnetic field, 
once in a direction parallel to the field, a second 
time in a direction perpendicular to the field. In 
the first case the electrical conductivity of the 
bar at liquid air temperature was 7.191073 
+0.02 ohm cm—, in the second case 6.18 10-3 
+0.02 ohm cm. After the bar had passed the 
transition temperature in the magnetic field the 
conductivity below this temperature was not al- 
tered perceptibly by the presence of the field, as 
could be established by changing its direction, 
switching it off, etc. The conductivity when the 
bar was cooled in the absence of a magnetic field 
was 6.45 X 10-* ohm cm, 

Hence there is a definite effect of the magnetic 
field orientation upon the electronic conductivity. 
The phenomenon can be explained by assuming 
that at the transition temperature the disorder- 
order transformation occurs preferentially in such 
a way that the tetragonal axis is that cubical axis 
of the high temperature crystal which lies nearest 
to the direction of the magnetic field. The 
electrical conductivity will then be highest in 
directions where Fe?+ and Fe* alternate, and 
lowest in directions where one meets either only 
Fe?+ or only Fe** ions. Hence it is to be expected 
that the conductivity will be highest parallel to 
the magnetic field ; that is, more or less parallel to 
the directions in which the c-axis has been 
“frozen in.” The picture of the disorder-order 
transition is in qualitative agreement with these 
experiments, which therefore corroborate the 
picture of Fe;O, given here. They show at the 
same time that the ionic picture according to 
which the electrons are thought to be added 
entirely to one-half of the Fe*+ ions is much too 
simple. A low transition temperature is very rare 
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for a crystal lattice in which the atoms are held 
together by strong bonds, and has generally only 
been found for lattices in which molecules are 
held together by weak van der Waals-London 
bonds. In the present case the transition can 
occur at such a low temperature because it is an 
electronic transition between different atoms. 
The low transition temperature also proves that 
the energy of order is very small. The entropy of 
order is R In 2 per atom, and accordingly TAS is 
of the order of 0.01 ev (T is of the order of 100°K). 
Neglecting the polarization and the displacement 
of the oxygen ions, one would calculate an energy 
of the order of a few electron volts (each Fe?+ ion 
is surrounded by 4 Fe** and 2 Fe?* ions in a first 
sphere of radius r=jav2, by 4 Fe*+ and 8 Fe*+ 
ions in a second sphere of radius r=}a/6, by 
6 Fe*+ ions at a distance 3av?2, etc.). The actual 
energy will probably be several (say 5) times 
smaller due to the polarization of the surrounding 
medium. This, however, is obviously far too high, 
and suggests that the wave. functions of the 
electrons extend over several atomic distances. 


2. MORE COMPLICATED SPINELS 
(a) Fe,AlO, 


The x-ray investigation! has revealed that this 
substance can be described as consisting for the 
main part of Fe(FeAl)O,, with Al in the octa- 
hedral holes. The intensity of the reflection with 
the lowest diffraction angle and the general rules 
derived in the preceding paper suggest that the 
Fe in the tetrahedral hole is the Fe** ion: 


Fe**(Fe?+Al)O, 


so that the substance belongs to the ‘‘inversed’”’ 
spinels, with an averaged distribution of divalent 
and trivalent ions about the 16-fold positions. 
Hence the substance appears to be analogous to 
Fe;0, 

Fe*+(Fe?+Fe*+)O, 


with the Fe** in the 16-fold position replaced 
by Al. , 
As already indicated by the x-ray intensities, 
the cation arrangement is not entirely that 
corresponding to the formula given above, and 
the best agreement with x-ray data was obtained 
by assuming that a small fraction of the Al has 
replaced Fe in the eightfold position. The sub- 
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stance was prepared by heating an intimate 
mixture of iron and aluminum oxides in a suitable 
gas atmosphere, and it is therefore always possible 
that at the high temperature of preparation the 
energy difference between the two arrangements 
is insufficient to outweigh the gain of entropy 
accompanying a more uniform distribution of the 
Al about both positions. 

These conclusions are corroborated by the elec- 
trical properties of the substance. If the substance 
were pure Fe*+(Fe?+Al)O, (or Fe?+(Fe*+Al)Ox,) it 
would be a rather poor electronic conductor, as it 
does not contain Fe*+ and Fe?* in equivalent 
lattice points. There might be a possibility of 
electronic transfer between Fe in the tetrahedral 
holes and Fe in the octahedral interstices, but 
this would only occur if the electronic levels in the 
two different positions were practically identical, 
which would be rather accidental. 

The electrical conductivity of the substance 
appears to be about 10* times smaller than that of 
Fe;0,: viz. 0.2 ohm cm-—. This is, however, sev- 
eral powers of 10 greater than would be expected 
for the pure substance, even if we take into ac- 
count that the sample showed a small deviation 
from stoichiometric proportions (for 1 Al ion it 
contained 1.005 Fe*+ and 0.995 Fe?+ ions). If the 
assumption that this residual conductivity is 
caused by a small fraction of Fe*+ in the 16-fold 
position is correct, one would expect the con- 
ductivity to be very sensitive to the heat treat- 


ment. Actually we found that by annealing the 
sample in a small (evacuated and closed) quartz 
tube at a temperature of 850°C for 4 hour, the 
conductivity decreased by a factor of 55, the new 
specific conductivity being 0.0037 ohm cm—. 
Obviously, as a consequence of this heat treat- 
ment at a temperature well below the formation 
temperature (1300°C), more Al has gone over to 
the position of lower energy, and the formula 
Fe*+(Fe?+Al)O, of the pure material has been 
more closely realized. 

This behavior is rather exceptional, as the 
conductive properties of other substances (to be 
described below) of similar type proved to be 
very insensitive to similar heat treatment, and 
such substances have even found technical appli- 
cation in cases for which a high thermal stability 
is required. 


(b) Solid Solutions of Fe;0, and MgCr.0, 


MgCr.2O, isa very poor conductor and can even 
be called an insulator (¢<10-® ohm cm). 
Solid solutions with Fe;0, were prepared in the 
following way. 

An intimate mixture of Fe2O3, MgO, and Cr20; 
in the proper proportions was pressed into bars of 
about 3 cm length and 0.3 to 0.4 cm diameter. 
These bars were preheated at 1200°C in a stream 
of pure COs, with the result that most of the Fe 
was reduced to the oxidation level of FesO4. Then 
the COz was replaced by very pure Ne and the 
temperature was raised to 1300°C. This tempera- 
ture was maintained for one hour. After that the 
carrier with the bars was shifted to a cooler part 
of the furnace and the substance was allowed to 
cool in the same atmosphere. For preparations 
containing more than 50 percent MgCr2O, the 
heating was conducted at about 1400°C, in order 
to obtain well-sintered bars. In this case the 
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atmosphere must be less reducing in order to 
prevent a deviation from stoichiometric pro- 
portions in the direction of FeO. Satisfactory 
results were obtained by heating at 1400°C in Ne 
containing 1 percent of Os. The content of ferrous 
ions was checked by dissolving a sample in 
chlorine-free hydrochloric acid under exclusion of 
atmospheric oxygen and determining the amount 
of Fe?* in the solution by a titration with 0.01 
normal ceric sulphate solution. 

The electrical resistance of a series of bars with 
varying composition was measured over a wide 
range of temperatures. The results are given in 
Table I and summarized in Fig. 2. 

These results show initially that there is a 
gradual transition from Fe;0., with a high 
conductivity and a low temperature coefficient, 
to MgCr2O,, with a very low conductivity and a 
high temperature coefficient. The curves of loge 
against 1/T are straight lines to a first approxi- 
mation, with a slope continually increasing from 
FesO4 to MgCr2O,4. The conductivity as a func- 
tion of temperature can be described for each 
sample by the equation 


o=0., exp(—e/kT), 


and the gradual decrease of the conductivity at 
room temperatures with increasing content of 
MgCr.O, is obviously mainly determined by the 
gradual increase of the “‘activation energy’”’ e¢, 
although the values of o.. found by extrapolation, 
formally the conductivities at T=, vary 
appreciably. 

The theory for the electronic conductivity of 
the oxides of these “transition elements”’ is still 
not very satisfactory.® It has been suggested? that 
the high conductivity of Fe;0, must be explained 
by the fact that, in an averaged distribution of 
Fe?+ and Fe*+ about equivalent lattice points 
(here those of the 16-fold spinel position), an 
interchange of electrons must be possible without 
requiring an appreciable amount of added energy. 
Hence the electrons would travel with a very low 
activation energy along the lattice points con- 
taining both Fe?* and Fe**, and the effect of the 
different ions Mg** and Cr** in solid solutions of 
Fe;0,4 with MgCr2O, must, as a first result, be 
_ 6]. H.de Boer and E. J. W. Verwey, Proc. Phys. Soc. 49, 
59 (1937). N. F. Mott and R. W. Gurney, Electronic 


Processes in Ionic Crystals (Oxford University Press, New 
York, 1940), p. 167. 
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the blocking up of this irregular Fe?+—Fe*+ 
pattern. It therefore follows that with increasing 
MgCr.O, content the electrons are forced to 
follow narrow and complicated channels to an 
increasing extent. These blocking barriers will be 
different from point to point in the crystal lattice 
of the solid solutions. Nevertheless, this dilution 
of Fe?+ and Fe** by the different ions of MgCr2O, 
can be described up to a temperature of about 
400°C with the aid of asingle “activation energy.” 
This behavior shows again that the conduction 
mechanism cannot be described completely with 
the aid of wave functions of the conduction 
electrons belonging to a single Fe*+ ion. On the 
other hand, the picture of the lattice wave 
functions (Wilson’s theory of semi-conductors) 
fails entirely to explain the absence of conduction 
in these oxidic substances containing the transi- 
tion elements Cr, Mn, Fe, Co, etc., and it is 
rather obvious that a better understanding of 
this case can be obtained only by introducing 
into the theory the interaction between a traveling 
electron and the positive hole left behind ; i.e., by 
going again in the direction of localized wave 
functions. 

Upon examining the curves of Fig. 2 more 
closely we observe a peculiar anomaly. It. is 
interesting to consider this point in some detail, 
as it is directly connected with the cation 
arrangement in these solid solutions and may 
throw some more light on the conduction 
mechanism in these materials. 

This anomaly is most easily demonstrated by 
plotting the activation energy of the different 
samples as a function of the composition (see 
Fig. 3). 

Starting from the building principles derived 
in reference 1 for the cation distribution in 
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spinels, we would expect’ that the different 
cations in these solid solutions would have the 
distributions about the cation lattice points 
shown in Table II (we only give the compositions 
with 25, 50, and 75 molecular percent MgCr20,). 

Hence from 0 to 50 percent content of MgCr20,4 
there are both Fe** and Fe?" ions in the 16-fold 
positions, and the electronic conduction will 
occur along planes containing both of these ions. 
At 50 percent MgCr2Ou,, however, the preference 
of Fe** for the eightfold position results in a 
separation of the Fe** and the Fe** ions. Hence a 
rapid increase of the “activation energy” can be 
expected if we approach this composition. There 
are a number of possibilities for the mechanism 
of electronic conduction for compositions con- 
taining more’ than 50 percent of MgCr2O,. This 
depends on the distribution of the divalent ions 
in the region 50 to 100 percent MgCr2O,, for 
which two possibilities are given in Table II, but 
for which an intermediate distribution between 
these two extremes might also be possible. X-ray 


TaBLE II. Arrangement of cations in solid solutions of 
Fe;0.+ MgCr.O,. 








Composition Eightfold position 16-fold position 


FesOu: MgCr2Ou (tetrahedral interstices) 


I 75:25 Fe*+ 
II 50:50 Fe*+ 
III 25:75 (a) }Fe**, 3Mg 
(b) Fe**, iFe**, {Mg 


(octahedral interstices) 


3Fe**, {Fe?*, {Mg, Cr 
3Fe?+, Mg, Cr 
iFe**, Mg, $Cr 

3Mg, $Cr 





evidence supports arrangement IIIa, but the 
differences in the calculated intensities corre- 
sponding to IIIa and IIIb are rather small, so 
that an intermediate distribution cannot be ruled 
out and is even highly probable in view of the 
small energy differences which can be expected to 
exist between the two extreme arrangements. 
Hence for compositions containing more than 50 
percent of MgCr2O, the conduction mechanism 
will be either an interchange of electrons between 
Fe*+ and a small amount of Fe** along the lattice 
points of the eightfold position, or an exchange of 
electrons between Fe*+ and Fe?*+ from different 
positions. In both cases the “activation energy” 
will be increased ; in the first case because of the 
blocking action of the Mg? ions and the large 
distance between the tetrahedral interstices 
(4av3 =0.43a), in the second case because of the 
difference between the energy levels at the 
octahedral and the tetrahedral lattice positions, 
respectively, although the distance is slightly 
shorter (4a11!=0.41a, both distances to be com- 
pared with the shortest distance between two 
octahedral interstices, }aV2=0.35a). The latter 
seems to be the less probable case. For example, 
in the case of AlFe2O, (Section 2(a)) there were 
no indications in favor of such a mechanism, 
which would only be possible if the difference 
between the two energy levels happens to be 
small. 


(c) Solid Solutions of Fe;0, and ZnCr.0, 


It is interesting to compare the behavior of 
these solid solutions of Fes04 and MgCr2O, with 
those of Fe30, and ZnCr2O,. Here the different 
ions will be distributed over the different posi- 
tions indicated in Table III. The Zn ions always 
occupy the eightfold position and are replaced by 
Fe*+ when the Fe;0, content is increased. 
Simultaneously the Cr ions in the 16-fold position 
are replaced by }Fe?++}Fe**. Accordingly, there 


TABLE III. Arrangement of cations in solid solution of 
Fe;04+ZnCr20.. 








Composition Eightfold position 16-fold position 


Fez:04:ZnCr2O.4 (tetrahedral interstices) (octahedral interstices) 





75:25 2Fe3t, 4Zn #Fe*t, 3Fe*+, 3Cr 
50:50 1Fe3+, 4Zn 4Fe2+, 4Fe3+, Cr 
25:75 1Fe3+, 4Zn iFe*t, 1Fe3+, 2Cr 
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are always equivalent amounts of Fe?+ and Fe*+ 
in the octahedral position, and there is no reason 
for a discontinuous increase of € with increasing 
ZnCr2O, content. 

The corresponding loge vs. 1/T curves are 
given in Fig. 4. At high temperatures the prepa- 
rations containing large amounts of ZnCre2O, 
show a break in the curves, suggesting that at 
these temperatures a second conduction mecha- 
nism covers the Fe*+ — Fe** conduction prevailing 
for lower temperatures. The curves can be 
represented by 


o=0., exp(—e/kT) +0.’ exp(—e’/kT). 


As the specific conductivity of ZnCr2O,4 is much 


greater than that of MgCr.20O, this phenomenon is 
most probably caused by the conductivity of 
ZnCr2O, itself, so that at higher temperatures the 
conduction is taken over by Zn (or by Zn and Cr 
ions jointly) as a consequence of the higher tem- 
perature coefficient of this conduction mechanism, 

In Fig. 3, therefore, where ¢ has again been 
plotted against. the composition, the value of e 
corresponding to the low temperature part of the 
curves has been used. The plot shows a continu- 
ous increase of e¢ with increasing ZnCr2O, 
content. 

We wish to acknowledge the valuable assistance 
of Mr. A. Bol in the preparation of specimens and 
the measurements. 
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Integral equations are derived that relate variations in the potentials of average force 


NE may define ‘distribution functions’’ for 

a given number, 7, of molecules as being 
proportional to the probability density of finding 
n molecules in a certain configurational position, 


between molecules of a system at two different densities or activities. These permit the cal- 
culation of the change in thermodynamic properties, or of the change in the distribution of 
molecules in space, in a liquid or crystalline phase, if either the temperature is varied, or if the 
the mutual potentials between the molecules is assumed to change. The equations are in a 
somewhat complex, but still distinguishable, matrix form. A matrix operates on the variations 
in potential occurring at one activity and transforms them into those occurring at a second 
activity. The matrix elements are combinations of the distribution functions at the second 
activity, to which the transformation is made, multiplied by powers of the difference of the 
two activities. The matrix appreaches the unit matrix in value as this activity difference ap- 
proaches zero. The product of the two matrices, one which transforms from activity @ to 
activity 8, with that which transforms from activity 8 to a, is the unit matrix. This leads to an 
integral (matrix) equation between the distribution functions at any two activities. 

The calculation of one element of this matrix product leads to a cell type equation for the 
computation of the activity, or free energy, of a system in terms of the potential between a 
single molecule and the molecules that immediately surround it in the system. 

The transformation matrices have solutions, at certain values of the activity difference y, 
corresponding to variations in the potentials of average force, which are transformed into zero 
variations at the new activity. The activities at which such variations can occur are those of 
the phase transitions in the systems. The solutions at these special values of the activity are 
differences of combinations of the distribution functions in the two pure phases that can coexist 
at these activities. 


INTRODUCTION 





averaged over all positions of the other molecules 
of the system. These functions are the exponent 
of minus the potential of average force between 
the molecules, divided by kT. The distribution 
functions of a given system at constant tem- 
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perature, but at two different densities, are then 
related to each other by equations which can be 
obtained by a simple extension of the concept of 
the Gibbs grand canonical ensemble. The special 
case of most general interest is that for which 
one of the two densities is the limit of zero 
density. At zero density the distribution func- 
tion for a classical system is the exponent of 
minus the potential energy divided by kT, and 
may be regarded as known. 

Although the equations relating the distribu- 
tion functions of a liquid or solid to those at 
infinite dilution are rigorous and exact, they are 
not susceptible of numerical solution, since they 
involve integration over the coordinates of the 
number of molecules in a macroscopic portion of 
the system. It is, however, possible to consider a 
variational method which leads to certain integral 
equations involving integrations over fewer 
numbers of molecules. 

The variational method discussed in this paper 
is essentially as follows: We consider the poten- 
tials of average force at one density, say those 
at infinite dilution, to be capable of some con- 
tinuous variation. This, of course, causes a vari- 
ation in the potentials at some other density, say 
that of the liquid. By assuming the variation in 
the potentials of zero density to be of the sort 
induced by a change in the mutual potentials 
between pairs, or between some relatively small 
number of molecules, the equations for the 
variation in the potentials of average force in the 
liquid can first be integrated over all but the 
molecules for which the change in potential is 
being computed, and the few with which these 
are interacting. One thus obtains a set of integral 
equations giving the variations in the liquid 
potentials as the sum of integrals of kernels 
times the variations in the potentials of zero 
density. The kernels are combinations of the 
unvaried distribution functions in the liquid. 

One of the purposes of this paper is merely the 
formal attainment of these equations. Their use 
is not further discussed in the body of this 
article, but will be dealt with in a subsequent 
paper. By considering variations due to changes 
in temperature, or other physically meaningful 
changes, one has a tool for the calculation of the 
thermodynamic properties of the system in terms 
of the distribution functions. 
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However, certain formal mathematical rela- 
tionships that follow from these equations appear 
to be of interest. The discussion of these rela- 
tionships occupies the main body of the paper. 


NOTATION 


Owing to the relatively complex equations 
involved, it is desirable to use a special symbolism 
in order to attain some brevity of expression. 
The final equations are in no way more com- 
plicated if a system of many components is con- 
sidered than if a single component system is 
handled. We deal with an infinite system con- 
taining o types of molecular species, and use the 
subscript s (or r) to refer to a special species of 
type s (or r). Symbols p, z, without subscripts 
refer to a set of quantities: 


*, Psy *** Poy (1) 
*. Re, ** *Re. (2) 


The symbol p, indicates the average concen- 
tration of molecules of species s, in molecules per 
unit volume, and 2g, is the activity or fugacity, 
defined by the equation: 


Zs= lim [pos exp(us—pos)/kT J, (3) 
p00 


P= Pi; Pa, °” 


2=21, 22, °° 


in which yp, is the chemical potential of species s 
in the system, and the limit is taken as all con- 
centrations, po, approach zero, with jo; the 
chemical potential, and po; the concentration, of 
species s, in this limit at which the system ap- 
proaches a perfect gas in properties. Thus 2, is 
proportional to the exponent of the chemical 
potential u,, divided by kT, with the propor- 
tionality constant so chosen that 2,/p, becomes 
unity in the perfect gas at infinite dilution of the 


system: 
lim p—0[z;/ ps | =1, (4) 


The letters k, m, n, x, pw, v, A, etc., when 
written without subscripts, will be used to indi- 
cate sets of integers: 


k=khk, ko, +++, Rs, ++ *Re, (5) 
of k, molecules of type s. The shorthand notation: 
R!i=k,!ko!---k,!-+-Re!, (6) 
k — ki nok2. ++ phe. sea ko 
p Pi “°p2 Os Po °; (7) 


zk =2)"igok- ° *Zah8+ ° * Ze", 


















































2m a ea nega 








se 
ex 
di 


es $ 
‘on- 
the 
1, of 
ap- 
Z, iS 
ical 
por- 
mes 
the 


(4) 


INTEGRAL EQUATION METHOD 189 


will be employed. With this notation the equa- 
tions for the « component system become iden- 
tical with those for a single component system of 
activity 2, or concentration p, in which k, m, n, 
etc., is used to indicate the number of molecules 
of the single component. The reader will probably 
find the ensuing discussion simpler if the one 
component case only is kept in mind. 

The coordinates of the ith molecule of type s 
are indicated by 


(ts) =Xie, Visy Zisy Qiisy si *Ofisy (8) 


including any internal coordinates, g, and the 
volume element, calculated in the Cartesian 
space of the constituent atoms by 


d(i,) —dxdyudsal “bag i a 
x 


with | the Jacobian transformation to the Car- 


tesian volume. The introduction of the internal 
coordinates causes no change in the ensuing 
equations. The reader will probably find the 
visualization of the operations easier if the case 
that no internal coordinates are present is kept 
in mind. The coordinates of k, molecules of type 
s are indicated by 


{ke} =(1e), (2.), i (Rs), (10) 
and the volume element by 
d{k.} =d(1,)d(2.) cae -d(ks), (11) 


whereas the coordinates of the set & of molecules 
are, EF 


{k} = {Ri}, {Re}, - {Re}, -+ {ke}, (12) 
and the volume element is 
d{k} =d{kijd{ks}---d{ke}---d{ke}. (13) 


A function, ¥,{”}, of the coordinates of the 
set {x} of molecules may happen to be simply 
expressible as a sum of functions of the coor- 
dinates {v}, of small subsets of molecules out of 





the set mw. The typical component function, 
¥,{v}n, composing a member of this sum, say for 
example, W2((7.)(js)), of the coordinates of the 
7,th and j,th molecule, will, if Y, is symmetrical 
in permutations of identical molecules, have the 
same functional dependence on the coordinates 
that the function W2((j,)(&s)) has upon the coor- 
dinates of the j,th and k,th molecule, but the two 
functions are distinct in the sum and have, in 
general, different values for some arbitrary value 
{2} of the coordinates of the set of molecules. 
The summation symbol, (>> {v},), will be used 
to indicate summation of functions of such sub- 
sets, {v},, over all subsets, for all values of », 
O<»v,<n,, including a single (constant) term for 
which all v,=0. It will be noted that there are 
n!/(n—v)!v! different distinct subsets, {v},, for 
each value of the set of numbers, v=71, ---, 
vs, ***, SO that in the trivial case that the func- 
tion is a constant, (—1)’, for all vs,*then 


(i {v}n)(—)’=0 if n>0, (14) 
({v}n)(—)’=1 if n=0. (15) 


Thus if Y, is one of a set of functions of the coor- 
dinates {nm}, symmetrical in permutations of 
identical molecules, defined for all values of 
including »=0, for which the function, Yo, is a 
constant, 


v,=V,{n}, (16) 
then we may define functions y by 3 
w= (Li{n}>)(—) "Wn, (17) 
and, in view of (14) and (15) we have that, 
Vn=(Ltv} add, (18) 
since then on inserting (17) in (18) we obtain 
the identity: 
Wn = (Li r} a) (Li}»)(—) 
= (LK) n)(Lte—K} nn) 
=V,. 


THE DISTRIBUTION FUNCTION EQUATION 


With this notation we now define distribution functions F,(p, T, {n}), of the coordinates of a set 
n of molecules, by the statement that in a system of infinite extent (V>>n,/p, for all s) at the tem- 
perature 7’, and concentration p, the probability that there are an appropriate set m of molecules at 
the coordinates {n} within the volume element d{n} is equal to p"F,(p, T, {n})d{n}. It is to be 
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noted that with this definition the functions F, are so normalized that, 


tim| — “ff; Fain) |=1. (19) 


We arbitrarily include Fy) =1 in the set of functions so defined. 

In a previous paper! it has been shown that for two sets of concentrations pq and pg, for which 
the pressures are P, and Pg, and for which the fugacity sets are z, and 2g, respectively, the dis- 
tribution functions are related by the equation set: 


Lexp(PaV/RT) | pa/Za]"Fn( pa, s 


ones m n+m 
-5< = af fo fe rovnal 2 =| Fasm(pp, T, {n-+-m})d{m}. (20) 


This actually represents two sets of equations for each pair of concentrations pa and pg, since obviously 
the subscripts on the right and left can be interchanged. Thus if we define: 








Gy aePav ier = F, (pa, T, {n}), (21) 

G,*=ePavieT| * | Fp, T, {n}), (21’) 
LZ, 

Z=Za— 2p, (22) 


we may write (20) as, 


Gum Ef fj: [orsnat, (23) 
Gt= ES tf + [ Garnd{m}. (23’) 


THE VARIATIONAL METHOD 


or as, 





Without regard, at first, to the physical interpretation of the functions G, and G,* we shall now 
proceed to carry out some mathematical manipulations which lead us to a set of integral equations 
that are of considerable interest in interpreting the relationships between these two sets of functions. 
The only mathematical limitations which must be imposed upon the functions G, and G,* are to 
require that they be symmetrical in permutations of identical molecules, and of such a form that 
the series (23) and (23’) converge absolutely. These conditions are assuredly fulfilled for the physical 
case. The symmetry in permutations of the coordinates of identical molecules is obvious. In addition 
we have that G, and G,* are positive definite, from their definitions in terms of probability functions, 
and for repelling molecules G, and G,* approach zero in value for all values of the coordinates {7} 
within a finite volume V as n/V approaches infinity. 

We now assume that G,, and therefore through (23’) G,*, depend on a parameter y, defining: 


V,=d InG,/dy, (24) 
=d 1nG,*/dy. (24’) 


For instance, we may take y to be the temperature 7, in which case, from (21) and (24) one could 
relate the. ¥’s to the various temperature derivatives of the system. The function F,(0, T, {m}), for 
p=0, is exp—(U,/kT), with U, the potential energy, if the system is classical. One might, for 


1W. G. McMillan, Jr. and J. E. Mayer, J. Chem. Phys. 13, 276 (1945) and for component systems, J. E. Mayer, J. 
Chem. Phys. 10, 629 (1942). 
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instance, consider that U, is U,°—yU,“™ with U,“ some sort of perturbation potential. In this 
case, if G,* referred to p =0, the functions Y,* would be 1/kT times the perturbation potential for 
the set of molecules. The functions Y, could then be related to the derivatives of the system 
with respect to this perturbation. 

We may now define the functions y and y* by (17) and by, 


= (Lo {n}»)(—) "Wn", (25) 
=(Liv}n)y*, (26) 


so that (18), and, 


are valid. 

We now proceed to derive Eq. (29) by which the function y can be calculated as a sum of integrals 
of kernels multiplied by the functions y*. The kernels are given by Eq. (30) as combinations of the 
functions G. A similar reciprocal set of equations, (29’) and (30’), can be used to compute the func- 
tions y* by integration over kernels formed from the G*’s multiplied by the functions y. These 
equations are obtained by differentiation of both sides of (23) with respect to y, obtaining: 


- “— . m>0 a 1 J Geena 
=r — (Fh tule) ff; + [ VissGend im} 
= z “Ztlo ff; fed x= Sf [G ee ei fe B} 


= EL f fj J veaGonatut. (27) 


In this the second line follows from the first by the use of (24’) and (26) and by considering that 


(DE {r} nm) =(Do fv} n) (Xo {w} m) =(Liv}ntu}m), 


where the last notation means summation over all subsets {v}, and {u}m including vy =0, and p =0. 
For any given value of the set m one may arbitrarily call m—y=k, and, remembering that there are 
m!/u!k! different subsets {u}m for any given value of the set m and y, one may multiply by this and 
sum over uw and k rather than m and {u}m. This leads to the third line. The fourth line follows by 
using equation (23). 

We thus have: 


> oof i fv Satu, (28) 


n>0 pw! 


To this equation we apply the operation of equation (17) to both sides, and obtain, 


w= “oi ff; + [VoemK usaf (m}, 


in which the kernels, Ki, m, are defined by the equation: 


Getnim a. 


30 
yoy (30) 


Kvn, nm= (Li K}x)(—)*™ 


The application of the equivalent series of operations to Eq. (23’), instead of (23), leads to the 
equations: 


vi= = NEU Lf og [tera wd tm, (29 
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with, 


Kian, nm=(L{«}x)(—)* * 





Special values of the kernels follow from (30) and (30’) to be, 


Koo0=Koo0=1, 

Kom=Gmn/Go; Ko.om=Gn/Go, 
Kian no=Krtnn0=0 if k>0, 

Knn0=Ks,,0=1. 


AN OVERSIMPLIFICATION 


Equations (29) and (29’) apply for all values 
of k>0. We may choose either the set yy, of 
functions or the set y;,* as arbitrary, limited only 
by the condition that they be symmetrical in 
permutations of the coordinates of identical 
molecules. 

A special case of the integral equation (29) (or 
(29’) but not both) has been used by Kirkwood 
and Monroe? to determine distribution functions 
for liquids. The special case chosen by them may 
be constructed from our equations as follows. 
We use (23) with the asterisk, G,*, referring to 
infinite dilution of a one-component system, and 
G, referring to the liquid state. Differentiation is 
then performed with respect to the coordinate 
of a molecule at a particular position over which 
the integration is not extended, that is, one of 
the coordinates of one of the molecules, (z), of 
the set ” in (23). The functions V,*kT are then 
the true forces, and are zero unless ” contains 
the molecule (7). Théy are symmetrical in the 
coordinates of all identical molecules except the 
molecule i. The y,* are assumed to be zero 
unless n = 2. The functions V,kT are the average 
forces in the liquid, which are related, by a well- 
_ known relationship to the unknown distribution 
functions G,, in the liquid. Equation (29) then 
has unknown y, and kernel K, but a subsidiary 
relationship between them is known. It then 
represents a non-linear integral equation, which 
Kirkwood and Monroe solved approximately 
with certain assumptions. 

The form (29’) for their problem would have 
presented them with a linear (non-homogeneous) 


2 John G. Kirkwood and Elizabeth Monroe, J. Chem. 
Phys. 9, 514 (1941). 


* 

Gitntm (30’) 
kn 

(31a) 

(31b) 

(31c) 

(31d) 





integral equation with known kernel K* and 
known function ¥*, but unknown average force 
y. The possible advantages of this form have 
not been investigated. 

The standard form of a non-homogeneous 
linear integral equation of the Fredholm type 
may be written as 


V(x) =¥*(x)—2 J K(xy)¥*(y)dy, (32) 


in which y and the kernel K(x,y) are known 
functions, and y* is the unknown function to be 
determined from the equation. There exists a 
reciprocal kernel K,*(x,y) depending on the 
value of the parameter z, such that, independent 


of y, 
W(x) =¥(2)-2 f KA(xy)¥(y)dy.  (32’) 


It is, of course, quite unimportant whether the 
variables x and y are one dimensional or a set 
of coordinates in multidimensional space. 

The reciprocal kernel, K*, varies continuously 
with the parameter z between the value z=0 
(where K*=K) and the first egienvalue, z,, for 
which the homogeneous equation, 


V@)=s, [KeyWord (33) 
has a solution. For values of z between 2, and 2;, 
the second eigenvalue, K* again varies con- 
tinuously with z. In general the functions which 
K* approach as z approaches 2, from above and 
below are different. 

Now let us consider a special case of (29), 
namely that there is only one component kind of 
molecule, and that ¥,*, which we can choose 
arbitrarily, is identically zero except for k=1. 
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There are, then, only two non-zero terms in the 
sum on the right, that with x=1, m=0, and that 
with x=0, m=1. We consider the case k =1, for 
which the kernel K;,1,9 when x=1, m=0 is unity 
from (31d). The equation for this case is, 


Wli) =v) +2 J Kio@, a*(d), (34) 


identical in form to (32). Were we justified in 
assuming that y,=0 for all values of k except 
k=1 then (29’) for k=1 would be, 


Yi(i) =yr(i) —2 f Kio) (avid)! (34’) 


of the same form as (32’). That is were the assump- 


° 2K ° 
tion correct then Ky, 9,1 and Ky, 9,1 would be simple 
reciprocal kernels for the parameter 2. 

We may now examine more closely what the 


kernels Ky,0,1 and Sia are. Choosing p3=0, 
z3/pg=1, in Eqs. (21) and (21’) defining G and 
G*, so that z.=z from (22), and using (30) and 
(30’) defining the kernels we have, 

F@), (i) 1 
—_— — K,{y), 35 


F, } 5 ; 
Ki0(@, w-f§e 


Kio1((i), (j)) = 


- F(a (35’) 


where the asterisks refer to infinite dilution and 
the functions F are the distribution functions at 
concentration p. Now if there are no internal 
coordinates then F,*=1, so that, 


Kies K—1, (36) 


which is a function symmetrical in permutations 
of the coordinates of the two molecules, ap- 
proaching zero in value as the distance between 
them becomes large. 

At least in the gaseous and liquid phases we 
also have F;=1, so that for these phases, 


p 
Ki,01=-LF2—1], (36’) 
Zz 


which is also symmetrical and approaches zeros 
for large distances. 

The simplicity of the result obtained from this 
unjustified assumption of identifying ¥,=0 for 
k>1 when y¥,*=0 for k>1 is made particularly 


intriguing when we consider that the function 
(36’) will show discontinuities in form as 2 
changes through the eigenvalues of the kernel 
(36), so that these eigenvalues would be iden- 
tified with the values of the activity at which 
phase transitions occur. 

The actual situation is far more complicated 
due to the occurrence of the summations in the 
original Eqs. (29) and (29’). We have discussed 
this simple case at some length because of its 
suggestiveness, and because it appears to repre- 
sent a simplified caricature of the true case. 


THE KERNELS AS OPERATORS 


We shall show that Eq. (29) may be regarded 
as a generalized matrix equation, in which a 
matrix, of which combinations of the kernels, K, 
constitute the elements, operate on the set of 
functions y,* to transform them to the set yx. 
The equation (29’) is the reciprocal equation, 
in which the reciprocal matrix, of elements 
formed from the kernels K*, operate on the set 
of functions y;, transforming them to the set y;,*. 

These matrices, which transform from the 
activity 2g to the activity Z., are functions of the 
distribution functions at the activity z., and of 
the difference, 2 =2.—2 8, Eq. (22), in the activity. 
We can thus consider the matrices, A,(y), as 
depending on two parameters, z and y, the matrix 
consisting of elements composed of the dis- 
tribution functions at activity 2, and trans- 
forming the set of functions y,,* corresponding to 
the activity z—y into the set y at z. We then 
obtain the matrix Eqs. (46) between the matrices 
for different values of z and y. 

Let us use the symbol a({k}, {m}) for a func- 
tion of the coordinates {k} and {m}, defined by, 


a({k}, {m})=(Li iw} elu} m) 
s+ 
5({«}x, {0} a)-—_$__—K i. am (37) 
(m—p)! 

where the summation runs over all subsets, {«},, 
of the set {k}, and all subsets, {u}m, of {m}, but 
5({«}x, {u}m)=0 unless the two subsets contain 
the same numbers of the same types of mole- 
cules, and unless the coordinates of all corre- 
sponding molecules of the two sets are identical; 
and in the latter case the integral of 6({«}x, {u}m) 
about the position of identity of the coordinates 
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is unity. We similarly define: 


a*({k}, {m})=(Qo xh ar (H} m) 


5({}x, tn a _ Ki; k,m—p+ (37’) 
(m—p)! 


With these two quantities we may write Eqs. 
(29) and (29’) as, 


Ve= LD J acth}, {m})yn*d{m}, (38) 


m>=0 
roe» forces, (om) oad im}, (38’) 
m>0 


These equations are now in a more conventional 
form than that of (29) and (29’). If @ and * are 
the sets of functions: 


D=yYo, V1, 7° Wky psitd (39) 
=yo", yi", , Vi", ita (39’) 


and A and A®* are the matrices of elements 
a({k}, {m}) and a*({k}, {m}), respectively; then 
(38) and (38’) may be written as: 


6=Ab*, o*=A*S, (40) 
AA*=A*A=E, (41) 


With E the unit matrix. The matrix multi- 
plication of the two elements involves both sum- 
mation over the middle index m and integration 
over the coordinates {m} for each value of the 
index m. The matrix A is thus a generalized type 
of matrix with continuous and discrete ranges. 
Thus we see that the Eqs. (29) and (29’) do cor- 
respond to the equations fortwo reciprocal com- 
mutative linear operators whose matrix elements 
can be represented by the a and a* of Eqs. (37) 
and (37’). 

Written in terms of these matrix elements 
equation (41) takes the form: 


7 ffs foc. {m})a*({m}, {n})d{m} 
“Zz JJ: far tk}, {m})a({m}, {n})d{m} 


=d({k}{m}). (42) 


This equation could be used, in principle at 
least, as a means of determination of the func- 
tions G from known functions G*. 


The matrices A and A* which we have defined 
by means of their elements in Eqs. (37) and (37’) 
depend both on the activity set z at which the 
kernels are determined (z, for A and 2g for A* 
in Eqs. (21) and (21’)), and upon the difference 
of activity, =2,—2s, between the two systems. 
Let us change our nomenclature at this place, 
and write the matrices as A,(y) of element, 


a(y, {R}, (m})=(Le tx} ala}m) 


ll 
5({«} 4 {ue} »)—_—— 
(m—p)! 
in which the kernels are made up according to 
(30) from the distribution functions at the 


————Ki,«m—n(8), (43) 


activity set z. 


If y becomes identically zero there remain in 
(43) only the terms for which »=m, and from 
(31c) and (31d) these are zero unless k= x=y=m, 
and are then unity. We thus have, 


limLA,(y) ]=£, (44) 


is the unit operator. 

Now let #(z) represent the set of functions 
Wo(z), ¥i(z), ---, ¥e(z) giving the variations in the 
potentials at the activity z. The transformation 
equations, (40), then take the form: 


&(z) =A,(y)P(z—y). (45) 
If we operate on this with A,,,(x) we find: 


&(z+x) =A,,2(x)®(z) =A.42(x)A-(y)P(z—y) 
=A,4.(x+y)P(s—y), 
or 
A.(x+y) =Az42(x)AA(y). (46) 


Equation (41) is seen to be a special case of (46), 
namely that for which x= —y, in which, from 
(44), A,+2(0) =E, and, 


Azi2(x)A—x) =E. (46’) 


A SPECIAL EQUATION 


Equation (42), which represnts a set of equa- 
tions, one for each value of k and n, is com- 
paratively: awkward to handle for the case of 
general values of k and n. It can be seen that 
the equation is trivially satisfied if n=0, as we 
shall show in this section. If R=0 and »=1, one 
obtains an expression which leads to an equation 
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for the activity, 2,, of the single molecule repre- 
sented in the set n=1. 

We use the symbol 7, or 7, to represent a single 
molecule of type s, and (i,) or (j,) to represent 
its coordinates, but use the subscript s to indicate 
a kernel in which the corresponding set m is a 
set consisting of one molecule of type s, and 
m—s when the set is some general set, m, 
(m,>1) less one molecule of type s. If ” repre- 
sents a set consisting of a single molecule of 
type s then 2"=g, in our notation. Using (37) 
and (37’), and Eqs. (31a) to (31d) inclusive, for 
special values of the kernels, we may write: 


a(0, {m}) =(2"/m!)Ko,o,m 
= (2"/m!)Gm/Go, (47) 
a*({m}, 0) =a({m}, 0) = 6(m, 0), 
a*({m}, (i,)) = —2Krm, 0,.+8({m}, (i,)). 
Using these, the sum }> a(0, {m})a*({m}, 0) 
has only the single term Go/Go=1 when m=0, 


so that Eq. (42) for the case k=0, n=0, is 
trivially satisfied. Similarly, all sums 


dX a({k} {m})a*({m}, 0) 


contain no non-zero terms if k>0O, from the 
second of these equations, so that (42) for this 
case is also trivially satisfied. 

The sum >} a(0, {m})a*({m}, (z,) =0, leads to 
the aie 


"a - Eas S- fGen adm 
0 m 


4, * 0. (48) 
23-——- => . 
c 


0 


If, in this, one sets the special case that the 
functions with asterisks refer to infinite dilution, 
see Eqs. (21), (21’) and (22), then s=z,, and 


2"Gm/Go= p” Fn(p, ie {m}), 


is the probability density that the set m of 
molecules will occur at the position {m}. The 


functions x 0,2 are now, from (30’), 


Fosel0, a { {4} m+ (ts) }) 
ee sae gee 


where the distribution function at zero concen- 











INTEGRAL EQUATION METHOD 195 


tration is: 
F,*(0, T, {u}) =exp—(U,{u}/kT), (49) 


and U, is the potential energy of the » molecules, 
if the system is classical. 
Equation (48) now becomes, 


(ps/2s) Fs(p, (is)) 


x m! ts i - [ Fale, T , {m}) 


X fexpl(—1/kT)(Uy4e{ {u} + (és) } 
—U,{u})]}dtuj}, (50) 


as a general equation for any system, if U is 
understood to be defined by (49) when the system 
is non-classical, in which case it is not a potential 
energy. 

If Eq. (50) is divided by V, and integrated over 
the coordinates d(i,), due to the normalization 
condition (19), one obtains p,/z,, the inverse of 
the activity coefficient, on the left-hand side. 

Let us assume a more specialized case, and one 
that is particularly simple. We assume that there 
are no internal coordinates, or that they have 
been integrated out of the distribution functions 
at zero density so that the potential U, is a 
function of the centers of mass, averaged over 
the internal coordinates. We further consider a 
single component system, and one existing in the 
liquid state so that F, is unity. If, now, it is 
legitimate to express U, as a sum of terms, 
uj, each dependent only on the coordinates of 
the pair of molecules, 7 and k, then if, 


fis=Lexp —(uij/kT) ]—1, (51) 
we have that, 


j=m 


Knoo= I fe (52) 


The functions f;; approach zero in value if the 
distance between the molecules i and 7 become 
large, since their mutual potential is then zero. 
Equation (50) now takes the form, 


=o ff. [Fn (p, T, {m}) 


j=m 


x II fid{m}. (53) 


j=l 
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The distribution functions F,, become zero if 
too many molecules are crowded into too small 
a space. The product of f;;’s is zero unless all of 


MAYER 


This equation is essentially a cell type equation 
in that one computes the activity, and therefore 
the free energy, by considering only the inter- 





the m molecules are within interacting distance actions of one molecule with those others that 
of the molecule 7. As m becomes large there are can approach it at any one time. The equation is 
no positions of the coordinates for which the exact in the form (50). It is not essentially dif- 
integrand is non-zero, so the sum converges, ferent from that derived previously by the author 
and, indeed, terms with m>13 must be very’ with Montroll,* although it is cast in an entirely 
small. different appearing form. 





EQUATIONS BETWEEN THE KERNELS 


To use Eq. (42) in the case of general values of k and , and convert it to a direct summation of 
integrals of products of the kernels K and K* is comparatively difficult. One runs into rather 
awkward summations, which then simplify with considerable labor. The simplest form of these 
equations can be found more readily by another method. However, we wish to emphasize that the 


equations derived in this section are actually another form of (42). 
We use Eq. (23’), write it for a set k++” of numbers, sum over yp, and divide and multiply each term 
in the sum by G;+,, using, however, (23) for G;4, in the numerator: 


Gu* si : ra | I f Gitn+y { 
2, 
K Grsntu : 
sf 2 SS f one 
> Grin Ludo v! 


gete Oe ee 
e ie» ——{— yf f-- f Tk+-n+pAk+p-+ yldi ut. 
u>0r>0 vip Gru 


In this equation we write »-+y4=m, and sum over all sets of numbers m. Since there are m!/v!u! 
different subsets, {4} m, for every value of the set of numbers m and yu, we may divide by this and sum 


over all {u}m. One obtains, 
“ie =f i J Chon (Ela ome Saal, [m}, 


- a 
Gipn= Gout m+k, k, n' d{m}. (54) 
x, 


If instead we start with (23), the equation, 


Giin= —ff- + [ Grimm m+k, k, n d{m}, (54’) 
x, m! ‘ 


may be obtained by the corresponding manipulations. 

At first sight Eqs. (54) and (54’) appear to be only complications of the relatively simple forms 
(23) and (23’), but the forms (54) and (54’) have one enormous advantage, namely that the sum- 
mations converge far more rapidly in these equations than in (23) and (23’). This superiority of 
convergence will be discussed in greater detail later, and is merely remarked here in passing. 

















* . 
Gisna= 





or with Eq. (30), 





3 J. E. Mayer and E. Montroll, J. Chem. Phys. 9, 2 (1941), Eq. (70). 
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We now use (54) to obtain an equation for Ka o,m, using (30’) defining this kernel as: 


+ waa o> (—s)" mae 
Khon = (Zlah)(—) Gta GE = Ef foo [Lhd awKnsce nll). 


However from (30’) we may write, 


Keen™(Z{A}es)Kitnne 
{u}. is any subset of x, and from (30), 


| rw aa™ (> { ut} «)Km+u, 0, ns 


where p= 


so that, 
(i x}x)(—)**Kuw wl tx, ca (> { xhafw}etA} ew) (—)**Kay, B mK mts, 0, n 


i (dX {u}e{A} ep) (—)e-Ae-r (SS { K—p—d} r-r-u( — jo) Rate, B, - am 0, n 


ss (> {uw} x) Ke, B, alate. 0, ny 


where the last line follows from the fact that the summation in the brackets of the preceding ex- 
pression is zero unless \+u=h. One therefore has, 


Las E Ete) f fog: [KaenKmica adi}. 


If, instead, (54’) is used, one finds, 


m20 Eta f f-; + [Kink m+x,0, nd{m}. 


, : + . 
Now for fixed values of the coordinates {nm} the kernels Kx o,, and K;,o,, are functions of the 
coordinates {k}, symmetrical in permutations of the coordinates of identical molecules. It is seen 


that Eqs. (55) and (55’) are identical to Eqs. (29) and (29’) with ee replaces ¥,* and Kz,o,n 


replacing yx. That is, the functions Kz, 0,» and K; o,n for any set of numbers u and any value of the 
coordinates {m}, regarded as functions of the coordinates {k}, form a suitable set y, and the cor- 


(55’) 


Ky, ee (55) 


responding solutions y;,*, respectively, to Eqs. (29) and (29’) or (38) and (38’). 


THE CONVERGENCE OF THE EQUATIONS 


Whether or not the sum in Eggs. (29) or (29’) 
converge depend on the variation functions yy; 
and y,*. We shall assume that the functions y,* 
are of the type which one might actually wish to 
handle if one considers variations in the potential 
energy of the molecules, or variations in tem- 
perature. In particular we assume that y,* 
approaches zero if the coordinates of any pair of 
the molecules composing the set k are sufficiently 
distant in the three-dimensional space of their 
mass centers, and that this approach to zero is 
more rapid than r~?, so that, 


tim| of f=; fyrate) |=finite (56) 





We shall also assume that the functions yy, 
approach zero in value as kK, or at least that 
the integral of Eq. (56) approaches zero as kR> 
We shall show that the form of the Eq. (30) 
defining the kernels is such, that, if the functions 
¥yx:* obey the conditions above, then if the 
kernels are formed at an activity for which the 
system can only exist in one phase, the Eq. (29) 
will converge and the functions y¥; will also obey 
the above conditions. If, however, the kernels are 
formed at an activity for which two phases are 
stable no such convergence can be expected. 
For such values of the activity set z that the 
system exists in a single phase, we may expect 
the distribution functions to show certain 
characteristics. These functions may be written 
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F,=exp(—W,/kT) where W, is the potential 
of average force. We would expect W, in single 
phase systems to have no terms corresponding 
-to long range forces. That is, we expect that for 
two sets of molecules, k and m, if the distances, 
Ym, between all molecules of the set & and all of 


the set m are large, then Wis» = WitWm. This 


means that: 
lim [Ferm ]=FiFm. (57) 


rm, k—> © 
Let us assume that Eq. (57) applies, and that the 
difference Fiim—F.Fm approaches zero more 
rapidly than 7;,,,-*, an assumption that cor- 
responds to having all forces approach zero more 
rapidly than 7r~%, 

Now let us examine the kernels, 


an F, n+rm 
Keenan (2) (Zleha)(-=, (68) 
z Frutn 

which equation is equivalent to (30) when p and 
z refer to the concentration and activity set for 
the functions without asterisks, see Eq. (21). 
Suppose the distances, 7:,n4m, between all mole- 
cules of the set k, and all molecules of +m are 
large. Then from (57) 


Fyangm = FF aim Frun=F Fa, 
and 
Fosuted Pete — Paul Pa 


But the sum (>> {x},)(—)** is zero if R>0, Eq. 
(15). It follows that, 


lim [KrtnnmjJ=0, k>O0, (59) 
Tk, n+-m— © . : 
and the approach to zero is more rapid than 


Yk,n¢m according to our assumptions. 
Let us now look at the product: 


* 
Krys n, n, mW n+m (60) 
which composes the integrand in the general 


term of the sum (29). If  — approaches zero as 
any of the distances of the set +m becomes 
large, that is it obeys (56), then the integrand is 
zero unless all molecules of the set »+-m are close 
together. But if the set ~+m is far from the set 
k, then from (59) the kernel is zero, and the 
integrand is zero. Furthermore, the argument 
leading to Eq. (59) applies if any one of the 
molecules of the set & is far from all of the set 


n+m. It follows that, within our assumptions, 
the integrand (60) is zero unless all of the 
molecules of the complete set k+n-+m are 
reasonably close together (within range of 
molecular forces one to the other) in the three- 
dimensional space of their centers of mass. It 
also follows, since we have assumed conditions 
leading to a sufficiently rapid approach to zero, 
that the integral of (60) over the coordinates {m} 
leads to a finite result independent of the volume 
V over which the range of integration is extended. 

If we also assume that ¥,*—-0 as kR-&, then 
we have assured the convergence of (29) since 
only a finite number of terms have non-zero 
integrands, and these each converge upon inte- 
gration. But we have also shown that, within our 
assumptions, each integrand is zero in (29) unless 
all the molecules of the set k are close together, 
so the y;’s also obey the condition (56). 

Now at such values of the activity set, 2, that 
two phases can coexist in the system, Eq. (57) 
is not correct. For instance, in the absence of a 
gravitational field (which would make z a func- 
tion of position in the system) the distribution 
function for the system would be: 


p"Fn=Xpa"F,%+(1 —x)p"F,™, 


where pa,p» are the densities of the two phases, 
and F, and F, the distribution functions in 
the pure phases. The number x is the volume 
fraction of the system in phase a, and is variable 
at the given z. Functions constructed according 
to this equation do not obey (57) even if F, 
and F,,® do satisfy it singly. We may conclude 
that, in general, solutions to (29) are not ob- 
tainable at such values of the set z that two 
phases can coexist in the system at the activities 
at which the kernels are determined. However, 
satisfactory solutions should be obtainable as 2 
is varied to approach the value of the transition, 
but the solutions will, in general, be different as 
this value is approached from the side of one 
phase or the other. 


EIGENVALUES AS THE ACTIVITIES OF PHASE 
TRANSITION 


The first terms in Eqs. (29) and (29’) are those 
for which m=0. From (31ic) and (31d) we see 
that only the term with x= in those of m=0 
has a non-zero kernel, and for this term the 
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kernel is unity. We may consequently write: 


wovit D (Dl «}e) 
m>1M: 


xf fj [KrenvEendio}, (61) 


* (—s)* 
=wt E (Eke) 
m>1 ™!: 


fg fis oternat ton, 61’) 


as modified forms of (29) and (29’). The fact that 
the terms independent of z in the matrix element 
were 6({Rk}{m}) has already been noted in the 
discussion under Eq. (43), and Eq. (61) and (61’) 
merely repeat this. 

If we define a matrix B,(y) by the equation: 


B.(y) =A.{y) —E, (62) 
so that, from (44), 
lim[ B.(y) ]=0, (63) 
y70 


then the transformation equation (45) becomes: 
#(z) = #(2—y) + By) ®(z—). (64) 
Let us write this for z=y, 
(2) = 8(0) + B.(z) (0), (65) 
and the reciprocal equation, 
(0) = @(z) + Bo(—z)P(z). (65’) 


Suppose that there exists a value of the set g, 
namely z,, for which there is a solution, ®,, to 
the homogeneous equation obtained by setting 
(0) =0 in (65’), namely: 


&,+Bo(—2,)’, =0. (66) 


This would mean that there would exist a pos- 
sible variation of the potentials of average force 
at the activity z, which would not be reflected 
in a change of the distribution functions at zero 
density. Without regard to the physical inter- 
pretation of this let us first examine its implica- 
tion in our equations. With this 4, for #(z,) in 
(65’) we find #(0)=0. However, if we insert 
(0) =0 in (65) we cannot obtain ®(z,) =4,40 
if the matrix Eqs. (61) (or (38)) converge. 

Now at such values of the activity set z that 
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the system exists in a single phase only, one 
expects that the distribution functions be 
uniquely determined, and that no stable vari- 
ations would exist that did not require altera- 
tions in the potentials at zero density. One also 
expects that no long range average forces be 
present, and that, therefore, as discussed in the 
preceding section, the matrix Eq. (61) would 
converge. 

However, at values of the activity set for 
which two or more phases are present the situ- 
ation is entirely altered. The distribution func- 
tions are not uniquely determined by the activity, 
since the relative fractions of the different phases 
vary at the constant activity. There are, then, 
variations possible which do not require changes 
in the potentials at zero density. But also the 
existence of more than one phase introduces a 
pseudo long range average force which com- 
pletely spoils the convergence of the matrix Eq. 
(61), corresponding to (65). 

We may therefore identify the values of z for 
which the homogeneous Eq. (66) has a solution 
with the values of the activity set for which more 
than one phase can coexist in the system. 

We have previously! identified the occurrence 
of a phase transition with the divergence of a 
power series in the activity. That is, one could 
start with a system of any activity, say activity 
z, and compute the thermodynamic properties 
such as pressure, density, or the distribution 
functions, F,, at some other activity, say z+y, 
in terms of the same property at z and a power 
series in y. If y were sufficiently small the power 
series converges, but if, as y increases, the 
activity z+y approaches the activity, 2,, at 
which a phase transition occurs, the power series 
diverges. 

Now the kernel, reciprocal to that in a 
Fredholm type integral equation, can be cal- 
culated by-an iterative process as a power series 
in the characteristic parameter, y, of the equa- 
tion, provided the value of y is less than the first 
eigenvalue of the corresponding homogeneous 
equation. As y approaches this value, the series 
diverges. 

One may regard Eq. (61’) as an equation for 
unknown functions ¥;, and attempt to solve it by 
using ¥,*=y, under the integral as a first ap- 
proximation, and use the resulting equation as a 
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second approximation. By reiterating this process 
one obtains an equation involving integrations of 
products of the kernels K* multiplied by y*, and 
by comparison with (61) one obtains an equation 
for K as a sum of integrals of products of the 
kernels K*. The sum is a power series in z. These 
equations, so obtained, can be transformed into 
the equations formerly obtained for the thermo- 
dynamic properties and distribution functions 
as a power series in 2. 

But the iterative process can only lead to a 
convergent series for those values of the param- 
eter 2 which are less than the first eigenvalue of 
the corresponding homogeneous equation. The 
two methods of obtaining the phase transitions 
are thus formally equivalent. 


THE SOLUTIONS OF THE HOMOGENEOUS 
EQUATIONS 

It is possible to find the significance of the 
solutions of the homogeneous Eq. (66) at the 
values, 2,, of the activity coefficient for which 
two phases coexist. 

We know, from experience, that the solutions 
G,, to Eq. (23), for »>0, undergo discontinuous 
changes as z crosses the values at which phase 
transitions occur. It is possible to disentangle 
the mathematical reason for this, and to see that 
the value of m for the term of maximum value 
undergoes discontinuous transitions at these 
values of z. In general, then, there exist two 
values, G,™, and G,™, which the sum (23) ap- 
proaches as z approaches z,, depending on 
whether the approach to the activity of the 
phase transition occurs from the activities for 
which phase (a) or phase (0) are stable. At z=2, 
any value, G,=xG,%+(1—x)G,® for 0<x<1 
is a physically meaningful solution to (23), cor- 
responding to having the volume fraction x of 
the system in the volume V in the form of phase 
(1). We might expect to use this for G, in (24) 
and take derivatives with respect to x to deter- 
mine the functions VW, for which no change in G,* 
is required. This procedure bogs down in mathe- 
matical difficulties. 

Let us instead turn to equations (55’). We 
again interpret the functions with asterisks to 
refer to infinite dilution. As z approaches z, from 


: ; (a) 
the two sides there are two solutions, _ o,n, and 


neh n, for every m and n, that satisfy (55’) with 


the same values of K*. It follows that: 


g Paci fs fre 


X[Ketcea—Ketncalim}=0. (67) 


The functions, 





(a) () 
Wht) = Kz on—Ki, 0, my (68) 


and the sets; 


®,((")) =p"), 7 ie Peri), (69) 


for any {m} are therefore solutions to the 
homogeneous Eq. (66) .at the activity z, at which 
phases (1) and (2) can coexist in the system. 


CONCLUSION 


The essential results may be summarized as 
follows. One may define distribution functions, 
F,, as being proportional to the probability 
density that a set, 2, of molecules are found at 


the coordinate postion {mz}. The logarithms of. 


these functions are —W,{m}/kT, in which 
W,{n} is the potential of average force between 
the molecules at the position {n}. These poten- 
tials depend on the activity set z of the different 
molecular species present, and vary with z con- 
tinuously from the true potential energy, at zero 
activity, but with discontinuities at the phase 
transitions. 

The function, W,(z, {x}) may be represented 
as a sum of terms, each depending on the coor- 
dinates, {v},, of a subset v of molecules of the 
set n, summed over all subsets: 


Wr=(Lir} n)wr(z,{¥} n)- (70) 


At zero density, at least, one usually assumes 
that only the functions w, for y<2, are non zero, 
corresponding to having the forces between 
molecules the sum of those between pairs only. 

One considers a possible variation in the 
various functions, w;(z, {k})/RkT. A set, &(z), of 
functions, Yo(z), vilz, (z)), ee vi(2, {k}), ae 
may be defined by 


Welz, {k}) =dLwn(p,{k})/kT I], (71) 


where 65(w;/kT) is the variation in w,/kT. This 
set of functions, depending on the activity 
parameter z, represents a self-consistent set of 
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variations of w;,/kT, of which the set at one 
particular activity, most naturally that at zero 
activity, may be regarded as arbitrary. 

In order to compute the set ®(z), from a knowl- 
edge of the functions at some other activity, say 
@(s—y), a transformation matrix A,(y) may be 


employed which transforms #(z—y) into ®(z). 


The elements, a.(y, {Rk}, {m}), of this matrix, 
given by Eqs. (43), (30), and (21), consist of 
combinations of the distribution functions at the 
activity z, to which the transformation is made, 
multiplied by powers of y, the difference of 
activities. The limit of the transformation 
matrix, as the activity difference y becomes zero, 
is, of course, the unit matrix. The matrix multi- 
plication involves both summation over all sets 
of molecules, and integration over the coordinates 
of each set. 

The product of the two matrices A,(y) and 
A,_,(—y) is the unit matrix. This is a special 
case of the more general relation that 


A,(y)Az+2(x) =Az42(x+y). 


These matrix equations are equivalent to a set 
of sums of integral equations, (55) and (55’), re- 
lating the distribution functions at different ac- 
tivities of the system. The equations are conver- 
gent if the activities are such that only single 
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phases are present in the system, even if the 
activity pairs are those at which different phases 
are stable. 

One member of these equations, (53), permits 
a computation of the activity or free energy of 
the system in terms of the interactions of one 
molecule with its neighbors in the system. 

If z is an activity for which the system exists 
in a single phase only, and 2, an activity at which 
two phases, (a), and (6), can coexist, then there 
are function sets ®, characteristic of z, which 
satisfy the transformation equation 


A,(z—2,)%,=0. 


In particular if z=0, Ao(—z,)#,=0, Eqs. (66) 
and (62), show that there exist a set, ®,, of 
variations of the potential of average force in 
the two phase system which can occur without 
variations in the temperature or potentials 
between the molecules at infinite dilution. The 
transformation matrix, A,,(y), to this activity, 
2,, of the phase transition, has elements such 
that the summations and integrations diverge. 
The solutions, ©,, to the homogeneous equation, 
Ao(—zy)®,=0, are differences of certain com- 
binations of the distribution functions of the two 
phases (a) and (b) which are at equilibrium at the 
activity zy, Eqs. (68), (30), and (21). 


4 
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Several formulae have been derived expressing the contribution to the intensity of electron 
scattering of restricted rotation in the vapors of ethane, ethylene, biphenyl, and their deriva- 
tives. The probability distribution of the restricted rotator about its axis of rotation is ex- 
pressed in terms of the barrier height, Vo. It has, therefore, been possible to evaluate the 
effect of barriers of known magnitude on scattering patterns. The possibility of evaluating the 
potential barrier using electron diffraction is discussed. 





T is well known from statistical thermodynamic 
studies! that many molecules exhibit re- 
stricted internal rotation at normal temperatures. 
This type of motion is of interest in the field of 
electron diffraction by gases because of its effect 
on the interpretation of scattering patterns in 
structure determinations, and also because of the 
possibility of using these patterns for evaluating 
the magnitude of the restricting potentials. Two 
methods are commonly used for interpreting the 
patterns from molecules which may possess re- 
stricted internal rotation. In both cases theo- 
retical intensity curves are calculated and are 
compared to the experimentally observed pat- 
terns. These curves are based either on a rigid 
model (which obviously neglects the internal 
motion) or on a model undergoing small vibra- 
tions according to the theory developed by R. W. 
James? for x-ray diffraction by gas molecules. 
This theory arises from a normal coordinate 
treatment which is generally quite unwieldy for 
polyatomic molecules and so its prescribed effects 
are introduced into the analysis in an empirical 
way. Another approach which has been tried? is 
to treat the molecule as a rigid structure with an 
additional internal librational degree of freedom. 
In this way the concept of a potential barrier 
hindering free rotation may be directly introduced 
and its effect evaluated. This procedure is justified 
in molecules possessing barriers of the order of 
3 kcal., as for instance derivatives of ethane, 
since in such cases the librational motion is far 
1A discussion of the method and a compilation of the 
results may be found in G. Herzberg, Infrared and Raman 
Spectra (D. Van Nostrand Company, New York, 1945). 
2R. W. James, Physik. Zeits. 33, 737 (1932). 
3 (a) J. Y. Beach and K. S. Palmer, J. Chem. Phys. 6, 639 


(1938). (b) J. Y. Beach and A. Turkevich, J. Am. Chem. 
Soc. 61, 303 (1939). 


more important than the other vibrational de- 
grees of freedom, especially at higher tempera- 
tures. On the other hand, some derivatives of 
ethylene may have barriers as high as 40 kcal. In 
this case the effect of the librational motion 
would be negligible at ordinarily attained experi- 
mental temperatures. In addition this motion 
and the energy quantum associated with it would 
approach the same order of magnitude as that of 
the other vibrations and could not be readily 
considered separately. 

P. Debye‘ has investigated the effect of the 
various types of internal motion, namely, vibra- 
tion and free and hindered rotation, on electron 
scattering by gas molecules. His theory for 
hindered rotation concerns 1,2 di-substituted 
ethanes and does not consider those cases for 
which classical statistics may not satisfactorily 
apply. In this article a quantum-mechanical 
distribution function for the angular distribution 
of the restricted rotator about its equilibrium 
position will be developed, and the theory ex- 
tended to include the more completely substi- 
tuted ethanes and also the derivatives of biphenyl 
and ethylene. For molecules undergoing hindered 
rotation, an increase in temperature may cause 
profound changes in their electron scattering 
patterns long before there are any considerable 
number of them undergoing free rotation, that is 
in states whose energy is greater than the barrier. 
This variation with temperature may be used to 
evaluate the hindering potential. 

The theory to be developed will apply to the 
case in which thé molecules of the gas occupy 
mainly the lower lying levels and the amplitude 
of the hindered rotation is not too large. The 


4P. Debye, J. Chem. Phys. 9, 55 (1941). 


202 


ELECTRON SCATTERING BY GASES 


other extreme, the case of free rotation, has 
already been investigated.* ® 


DEVELOPMENT OF THEORY 


The expression for the intensity of coherent 
molecular scattering from a rigid molecule is 


Sif; sinsr;;/sriz, (tJ) (1) 


i=1 1 

where NV is the total number of atoms in the 
molecule, f; is the atomic scattering factor for 
electrons for the ith atom, 7;; is the distance be- 
tween the ith and jth atoms and s=47 sin(@/2)/d 
where @ is the angle of scattering and X is the 
wave-length of the electrons. We are concerned 
with the contribution from varying interatomic 
distances which may be expressed by integrals of 
the type 


f sin(s7;;)P(rij)drij/srij. (2) 


The interatomic distances, r;;, affected by the 
internal rotation and their probability distri- 
bution functions, P(r;;)dri;, will be expressed in 
terms of the angle of internal rotation, ¢. By 
solving the Schroedinger equation for the re- 
stricted rotator undergoing small vibrations, the 
probability distribution of the rotating group as 
a function of angle ¢ for each of the lower states 
can be found. The distribution of the gas 
molecules among these states is determined by 
the Boltzmann distribution law. The probability 
functions P’(¢)d¢=P(r;;)dr;; may be obtained 
from multiplying the probability distribution 
within each of the hindered rotational states, as 
obtained from the Schroedinger equation, by a 
factor governing the frequency of occurrence of 
each of these states, as obtained from the 
Boltzmann law, and adding these products. 

The effect of hindered rotation on electron 
scattering will first be developed for substituted 
ethanes of the type X3;CCY3. They are symmetric 
top molecules and therefore lend themselves 
readily to theoretical treatment. It has been 
shown® that for symmetric top molecules the 
Schroedinger equation for rotational motion is 
separable into an equation containing the vari- 


ed: Karle, J. Chem. Phys. 13, 155 (1945). 
°H. H. Nielsen, Phys. Rev. 40, 445 (1932). 
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ables governing the motion of the molecule as a 
whole, and one in terms of the variable concerned 
with the internal motion. We are interested in the 
latter which is 


Py /dg? +82? I(E— V) p/h? =0, (3) 


where J is the reduced moment of inertia equal to 
I,I2/(Ii1+J]e), where I; and I; are the moments of 
inertia about the axis of the molecule of the two 
ends rotating with respect to each other, and ¢ 
is the variable of internal rotation. The other 
quantities have the usual significance. The po- 
tential function, V, is assumed to have the form 


V = Vo(1—cose¢) /2, (4) 


where Vo is the barrier height and o is the number 
of minima in the potential function, for instance 
three for ethane. When ¢ is small, we get 


V=V(0¢/2)’. (5) 


On substituting (5) into (3) the Schroedinger 
equation for the simple harmonic oscillator is 
obtained. It is 


dy /dg? +82 ILE — Vo(o¢/2)? y/h?=0. (6) 
The normalized solutions to this equation are? 


¥n(o) = Nn exp(—ad’/2)H, (a4), (7) 
where 
Na={(a/m)*/2"!}3, 
=20°r IV o/h’, 


and H, is the Hermite polynomial of order n. 
Since we are concerned with the situation in 
which only the first few states are occupied 
significantly, we may add to the distribution 
function for these lower states the distributions 
for the higher harmonic oscillator states since 
their Boltzmann factors are so small that their 
contribution is essentially equal to zero. The 
advantage in including them arises from the fact 
that the infinite sum may be readily written in 
compact form. 

The probability distribution function in terms 


7 There are also vibrations of small amplitude which may 
be considered as centered not only about ¢=0, but also 
about ¢=22n/o¢ as n goes from 1 to o—1. These lead to 
wave functions and a distribution function which assume 
significant values in the vicinity of these « values of ¢. For 
our purposes the desired results may be obtained from 
considering only the distribution about ¢=0 since the 
others make an identical contribution to the scattering 
problem. 
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of the contributions from the various states is 
(a/m)*(1—expl —hv/kT ]) exp(—a¢’) 
X>¥ An?(ak) exp(—nhv/kT)do/2"n!, (8) 


n=0 
where v=o(V o/2J)/2z, T is the absolute temper- 
ature, and k is the gas constant. The quantity 
1—exp[l—hv/kT] normalizes the Boltzmann 
factors. Mehler®in studying the parabolic cylinder 
function, D, obtained the generating function, 


exp[(—[1—#¢]/[1+¢])p?/2]/(1—#)! 
pie "D,*(p)/n!. (9) 


n=0 
The following relation exists between the para- 
bolic cylinder function and the Hermite 
polynomials, 


D,7(p) =2-"Hy?(p/2!) exp(—p?/2). (10) 
On substituting (10) into (9) it is found that 


i) 


X PH,%(p/28) exp(—p*/2)/2en! 
=exp[—((1—¢]/[it+eDe%/2/(—e). (11) 


By comparing (11) to (8) it is apparent that by 
setting t=exp[ —hv/kT] and p/2'=a'¢ we get 
for the sum of (8) 


(a/r)*[tanh(hv/2kT) } 
Xexp[ —a¢? tanh(hv/2kT) |d¢. 


This is the normalized distribution func- 
tion we have been seeking. The function 
b exp[ —b*¢? |d¢/z', where 6? =a tanh(hv/2kT), is 
the probability of finding the rotating group 
between ¢ and ¢+d¢. Expression (12) is equiva- 
lent to the distribution function obtained by 
Bloch? for the linear harmonic oscillator. If hy is 
small with respect to RT we obtain the classical 
result from (12), 


(o2Vo/4ekT)* exp(—o2Vod?/4kT dd. (13) 


We may proceed now with the scattering 
problem and the next step is to express the 
interatomic distances which change with the 
internal rotation in terms of angle ¢. The equi- 
librium distances which occur depend upon 
whether the molecule, X;CCYs, is at equilibrium 


(12) 


8 F.G. Mehler, J. f. reine u. angew. Math. 66, 161 (1866). 
°F. Bloch, Zeits. f. Physik 74, 295 (1932). 
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in the staggered or eclipsed configuration. The 
staggered configuration is the more probable one 
and is considered first. Two types of distances 
appear, 

r1= (ci +c cos¢)! 
and (14) 
r2=(c1—C2 cos|¢—7/3 ])}, 


where ¢,= (l2?+1,;")/2, and ce=(l.2—1,?)/2. The 
quantity J, is the smallest possible distance of 
approach of atoms X and Y if the rotating groups 
could rotate freely and J, is the largest possible 
distance between them. 

The contribution from small variations in 7; 
resulting from hindered rotation is obtained from 
evaluating the integral, 


| sin[_s(ci+cs cos)? | 
Xexp(—b262)do/m!s(c1+c2 cosd)?. (15) 
We express cos¢ as 1 —¢?/2 and find that 


sin[s(ci+ce cos)? | 
~I.P. exp[i(sle—sced?/4l2)] (16) 


where I.P. refers to the imaginary part. Integral 
(15) becomes the imaginary part of 


b explisls) { exp[ — (b?+7B)¢? |do/m'sl2, (17) 


where B =sc2/4l2. The small variations in 7; have 
been omitted in the denominator of the integrand. 
Integral (17) has the value b exp(isl2) /sl2(b?+7B)! 
and the imaginary part is 


b[cos(sl2) I.P.(6?—7B)? 
4-sin(sls)R.P.(6?—iB)*]/slo(b4+B?)}, 


where R.P. refers to the real part. By setting 
(b?9—1B)*=p' exp(zé/2) it can be shown that 


I.P.(6?—iB)' = — (64+ B?)! sin(é/2) 


(18) 


and (19) 


R.P.(6?—7iB)' = (64+ B?)! cos(é/2), 
where tant = B/b? =sc2/4l.b?. We thus obtain for 
the contribution from 71 

(cosé)? sin(slp—£/2) /sle. 


There are three r; terms occurring for X3CCYs. 
If only a rigid structure was considered, the 
contribution to the scattered intensity would be 


(20) 
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sin(sl2) /sle. The effect of the restricted rotation is 
to shift the maxima and minima by angle £/2, 
which varies with s, and to introduce a damping 
factor, (cost)#. When the same atoms oscillate at 
each end of the C—C bond the relation ],?—1],? 
=4q? holds where a is the radius of rotation. 
When this relation is substituted into (20) a 
result similar to Debye’s is obtained. 

The contribution from small variations in re 
resulting from hindered rotation is obtained from 
the integral 


if sin[s(c1—c2 cos_¢@—7/3 ])*] 


Xexp(—b*¢?)do/m's(¢1—colo—2/3])3. (21) 


For small values of ¢ we have cos(¢—72/3) ~0.5 
+(3'p) /2—0.25¢?. If we define 


T1 =$(¢1—C2/2)}, 
B=C2/8(¢1—C2/2) —3(¢2/[¢1—c2/2 })?/32, 
and 


y= 33c2/4(c1 —C2/2) 
(21) becomes 


if sin[71(1—y¢+6¢?) ] 
ts Xexp(—b%$2)db/a4n. 


(22) 


As in (17) the variation of r2 about its equilibrium 
value is omitted from the denominator. Integral 
(22) may be expressed as 


inns f cos(y71¢ — Bri¢”) exp(—b*g?)do 


—oO 


—_ cosrs f sin (y71¢ — Bri?) 


—o 


xexp(—198)d6| /n, (23) 


After further expanding the integrands and 
bearing in mind that odd integrands do not 
contribute to the values of the integrals, it is 
found that the problem reduces to evaluating 


J costes) expl-('+i8ri)o*]d¢. (24) 


—e 


The real part of the integrand of (24) is the 
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contributing portion of the integrand of the first 
integral in (23), and the imaginary part is the 
contributing portion of the integrand of the 
second integral in (23). Integral (24) has the 
value 


m exp —y’r1"/4(0?+7671) ]/(0?+7671)*. (25) 
By defining (b?—787:)!=p' exp(—70/2) we .find 
that the real part of (25) is 
m cos[_(y?8r1*/4[b4+- 671? }) — 0/2] 
Xexp[ — 2b?71?/4(b4 +671") ]/(b4+- 671"), 
and the imaginary part of (25) is 


a sin[ (y?67r13/4[b++6?7 1? ]) —6/2 ] 
Xexpl —°b?ry?/4(b4+ B71’) ]/ (t+ fr")! (27) 


where tan@=£7;/b?. By combining (26) with the 
coefficient of the first integral in (23), and (27) 
with the coefficient of the second, we get the 
desired result, the contribution from re, 


(cos@)! sin[71— (y?71 sin?0/48) +0/2 ] 
Xexp[ —~’r1? cos?0/4b? ]/71. 


(26) 


(28) 


If distance rz were fixed at the equilibrium point, 
the contribution to the scattering would be 
sint;/71. The hindered rotation introduces not 
only a shifting in the maxima and minima and a 
damping term, (cos@)?, but also an exponential 
damping term dependent upon s?. There are six 
distances of the type 72 in XsCCYs. 

When equilibrium occurs in an X;CCYs; 
molecule in the eclipsed position, there are again 
two different types of interatomic distances 
affected by the hindered rotation, 


r3=(c1—C2 cos@)! 
and (29) 


%= (c1—co[ o— 2/3 ])?. 


The contribution to the scattered intensity may 
be obtained in precisely the same fashion as for 
the molecule in a staggered configuration at 
equilibrium. The contribution from 1; is 


(cosf)? sin(sli+¢/2)/sl 
where tan{ = sc2/4b*l,. The contribution from 7, is 


(cosé)? sin[r2+ (u?re sin?6/4w) —6/2 ] 
Xexp[ —u?r2? cos?6/4b? ]/72, 


(30) 


(31) 
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where 


tané = wr2/b’, 
w = €2/8(¢1+62/2)+3(¢2/[e1tc2/2 })?/32, 
72=S(¢1+¢2/2)}, 
and : 
w= 3'c2/4(c1+¢2/2). 


Molecules of biphenyl and its derivatives, 
ethylene and its derivatives, and partially substi- 
tuted ethanes also exhibit hindered rotation. In 
general those molecules are not symmetric tops 
and therefore the internal rotational motion 
cannot be separated from the over-all rotational 
motion to give Eq. (3). If the coupling is not too 
strong, and this may be determined by comparing 
the values of the principal moments of inertia, it 
is still legitimate to consider (12) as the correct 
distribution function arising from the lower lying 
levels. 

For many bipheny] derivatives the equilibrium 
position probably occurs with the planes of the 
benzene rings at right angles. In this case the 
formula to be used is 


sin(sc;') exp[[ — s¢o?/16c,b? ]/sc1}. (32) 


If the planar configuration is correct for some 
biphenyl derivatives as it is for ethylene deriva- 
tives, formulae (20) and (30) apply. It should be 
born in mind that the two potential minima may 
not be the same. 

In substituted ethanes the structure problem is 
rather complicated owing to the existence of a 
threefold potential which does not have identical 
minima.!® For the case of HaXCCXHag, the trans- 
configuration has been found to be the more 
stable one.* The restricted rotation formula 
which applies to this model is (20). It is probably 
more nearly correct in using (20) to set o=3 
instead of one, even though it does not have 
identical minima. The potential barrier Vp would 
still signify the energy required to break out of 
the trans-position, but its value might be different 
than that of the potential barrier at the cis- 
position. 

The formulae for hindered rotation occasion- 
ally may not be reliable for Vo/RT less than 5 but 
often may be used for Vo/RT as low as 3. This 
estimate is based on the requirement that the 


10 For an excellent discussion of this problem see reference 
1, p. 346. 
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Fic. 1. Variation of scattering pattern with changes in 
Vo/RT for Cl;CCCl; and HeCICCClHe. Vo/RT is equal 
to © in A and D, 5 in B and E, and zero in C and F. The 
hydrogen distances in H2CICCCIH2 were not included in 
the calculations. The abscissa is given in s units. 


distribution function P’(¢) should damp out fast 
enough to justify the replacement of 1—cose¢ by 
(o@)?/2 in (4), and also damp out fast enough to 
allow the approximated interatomic distances as 
in (16) and (22) to yield valid results when s 
becomes as large as 20. 

It is interesting to compare the effects of 
restricted rotation on distances whose equilibrium 
positions are trans, cis or intermediate. Trans- 
and cis- formulae are represented by (20) and (30), 
respectively, whereas (28), (31), and (32) repre- 
sent various intermediate positions. The ex- 
ponential damping factor in the latter three 
formulae should be noted. Owing to this term the 
contribution to the scattering from intermediate 
distances affected by hindered rotation is largely 
damped out, whereas that from the cis- and trans- 
positions is relatively unchanged. This may be 
related to the case of free rotation in which the 
contributions to the scattering arise from the 
vicinity of the cis- and trans-positions. 


APPLICATION OF THEORY 


Electron diffraction offers the possibility of 
determining the height of the potential barrier 
restricting rotation in certain molecules. One 
method, which is probably the least promising, 
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involves treating the barrier as an additional 
unknown in fitting theoretical intensity curves of 
assumed models to observed data. The difficulty 
with this method is that the number of variables 
ordinarily involved is too great to permit the 
establishment of a unique solution. Another 
method which is based on the radial distribution 
function has been suggested by Debye.‘ He has 
shown how the values of the radial distribution 
function in the vicinity of an equilibrium distance 
can be related to the probability distribution for 
that distance, e.g., (12). 

Another approach to the problem of evaluating 
the barrier, which obviates to a large extent the 
ambiguity which arises in the first method sug- 
gested, is to study the changes introduced into 
the diffraction patterns when the sample of vapor 
is heated to successively higher temperatures. In 
favorable cases the ambiguity may be resolved 
owing to the fact that only a unique combination 
of structural parameters and potential barrier can 
fit the diffraction data over a range of tem- 
peratures. 

Theoretical intensity curves for various 
values of Vo/RT are plotted for Cl;CCCl; and 
H.CICCCIH, in Fig. 1. Curves A, B, and C refer 
to Cls;CCCI; and were based on an assumed model 
in which C—C=1.54A, C—Cl=1.76A and angle 
CI—C—Cl was given the regular tetrahedral 
value. Curves D, E, and F refer to HoCICCCIH, 
and were based on the same model as hexa- 
chloroethane by giving angle Cl—C —H the same 
value as was given to angle Cl—C—Cl. Curves A 
and D were calculated from the scattering 
formula for a rigid molecule, (1), the atomic 
scattering factors being set equal to the corre- 
sponding atomic numbers. Curve A was based on 
a staggered configuration and curve D on a trans- 
configuration. These curves represent the scat- 


tered intensity when the barrier is very large or 
the temperature is very low. Curves B and E 
were calculated for Vo/RT equal to 5. The 
contribution from the rigid portion of the mole- 
cules was calculated using formula (1). The 
contribution from the hindered rotation was ob- 
tained using formulae (19) and (27) multiplied by 
the atomic number of chlorine squared, and the 
proper frequency of occurrence of these terms for 
the molecules considered. The quantum correc- 
tion was not introduced and o was set equal to 
three. Curves C and F represent free rotation in 
the molecules. Formula (1) was used to calculate 
the contribution from the rigid portions of the 
molecules and the contribution from the freely 
rotating portions was calculated from the for- 
mulae in reference (5). Additional correction 
terms to the asymptotic formula of Debye were 
found to be negligible. 

The effect of the internal motion on the scat- 
tering patterns of hexachloroethane and _ 1,2- 
dichloroethylene is vastly different. A great 
change in pattern has occurred in the range © to 
5 for Vo/RT in hexachloroethane whereas practi- 
cally nochange has occurred for dichloroethylene. 
If T is about room temperature then Vo/RT 
equal to 5 means that Vo is equal to 3 kcal./mole. 
This is a customary order of magnitude. 

We may conclude that electron scattering from 
some molecules may be sufficiently sensitive to 
changes in Vo/RT to permit an evaluation of 
their potential barriers. The possibility of a 
useful experiment depends upon the magnitude 
of the barrier, the vapor pressure and the heat 
stability of the substance as well as the pattern 
sensitivity. It is perhaps well to emphasize that 
the determination of the structure of some 
molecules may be greatly dependent upon proper 
evaluation of their hindered internal rotations. 





THE JOURNAL OF CHEMICAL PHYSICS 


Letters to the Editor 








HIS section will accept reports of new work, provided 
these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 





Effect of Temperature on the Mass Spectra 
of Hydrocarbons 


R. E. Fox anp J. A. H1PPLE 
Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 
February 26, 1947 


EVERAL years ago one of us reported a marked tem- 

perature effect on the relative peak heights observed 
in the mass spectra of hydrocarbons.! As this work has 
not previously been published, it would appear desirable 
to describe the experiment now in view of the recent 
interest expressed in it privately. 

Figure 1 shows the ion source of the mass spectrometer 
tube used in the investigation. The thermocouple on the 
gas inlet block (3) should provide a fairly good measure 
of the temperature of the gas emerging from the small slit 
in it and crossing the electron beam. The ions formed by 
electron impact are then accelerated through the slits in 
electrodes (2) and (1) and pass into the conventional 
analyzer of the sectored-field type. The temperature of (3) 
was varied by adjustment of a furnace wrapped around 
the glass envelope of the ion source. 

Using a mixture of isobutane and neon, the intensities 
of various ions in isobutane relative to Ne®+ were measured 
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Fic, 1. Arrangement of ion source for studying the effect of temperature 
on the mass spectra of hydrocarbons. 
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M/2g AS A FUNCTION OF TUBE TEMPERATURE 
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Fic. 2. Variation with temperature of several ions observed in isobutane 
(neon used as internal standard). 


for various values of the temperature of (3). These data 
are plotted in Fig. 2. It is seen that there is a very marked 
dependence on temperature of the C,Hio* ion. In fact, this 
ion is almost twice as abundant at 70°C as the C4Hg"* ion, 
becoming equal to it at approximately 180°C and is rela- 
tively less abundant at higher temperatures. Thus, the 
ratio of these two peaks may be used to measure the tem- 
perature of the source. The mass 58 peak in n-butane has a 
similar temperature dependence. Other hydrocarbons show 
an even greater variation of the spectrum with change of 
temperature. For instance the parent peak of 2,2,3-tri- 
methylpentane is more than five times as large at 175°C 
than it is at 225°C. Other peaks in this spectrum show 
changes that are less pronounced. 


1J. A. Hipple, reported at the St. Louis meeting of the American 
Institute of Chemical Engineers, November, 1944. 





Decomposition of Solid Barium Nitrate 
by Fast Electrons* 


AUGUSTINE O. ALLEN AND JOHN A. GHORMLEY 
Clinton Laboratories, Oak Ridge, Tennessee 
February 21, 1947 


HE behavior of irradiated salts on solution is an 

almost unexplored subject, which has become of 
interest recently in connection with production of radio- 
isotopes. 

During 1944, the authors had occasion to irradiate ordi- 
nary C.P. grade barium nitrate crystals with a beam of 
1.2-Mv electrons from a Van de Graaff generator. On dis- 
solving the irradiated salt in water, copious effervescence 
occurred, and nitrite was found in the solution. 

In one experiment, the salt was irradiated in a small 
closed glass vessel under helium; the energy dissipated was 
1.1 watt-hr./g of salt. After irradiation, the vessel was 
evacuated and water was admitted; gas was evolved on 
solution to the extent of 0.9 cc/g of salt. Analysis showed 
the gas to contain 3 percent of hydrogen; the balance con- 
tained no acid gases and was presumed to be oxygen. The 
solution contained nitrite ion. 

In a second experiment, the salt was irradiated while 
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TABLE I. Formation of gas and nitrite on solution of barium nitrate 
irradiated with 1.2 Mv electrons at a rate of approximately 1.5 watts/g. 








Gas formed, Calc. Oz equiv. 
Irrad. % NOs~ converted cc/g (22°C, to NOz~ found, 
time, hr. to NO.- approx.73cm) cc/g (22° C, 73 cm) 





$$ 4. 
10.5 , 9: 
13.0 - 1, 
18.0 ‘ 6 
21.1 8. 








spread out on a pan in the open air. The crystals initially 
turned yellow, but on further irradiation became milky, 
so that the color was obscured. Samples were taken from 
time to time and placed in a closed tube, to which water 
was admitted after evacuation. The gas was not analyzed, 
but the volumes produced were measured, and the nitrite 
content of the solutions was determined. Results of the 
measurements are shown in Table I. Energy was not 
absorbed uniformly by the salt during the irradiation 
because of variation of thickness of the salt layer and its 
distribution on the pan with respect to the electron beam; 
the rate of energy dissipation was estimated as roughly 
1.5 watts/g salt. Solutions of irradiated salt (0.01-0.07M) 
showed a pH of 9 to 10, and were thus more alkaline than 
can be accounted for by hydrolysis of nitrite. 

At 25 watt-hr./g, about 21 percent of the nitrate had 
been converted to nitrite. The conversion rate decreased 
with time, suggesting that a steady state might eventually 
be reached after 40 percent or more conversion. 

Qualitative experiments showed that gas evolution on 
solution was not greatly reduced by heating the irradiated 
salt below the melting point, but was entirely prevented 
by fusion of the irradiated salt. Presumably the gas came 
out during fusion, but no direct observation was made on 
this point. Irradiation lowered the melting pcint of the 
salt, presumably because of the admixture of nitrite. 

Narayanswamy! has reported the formation of nitrite 
near the surface of nitrate crystals illuminated by light of 
wave-length 2200-2500A. He ascribed the localization of 
the reaction at the surface to the escape of oxygen gas, 
assuming that recombination of the products should occur 
in the interior. In the present experiments, practically all 
the gas remained trapped in the crystal (as O atoms or O» 
molecules), so that thermal recombination must occur 
very slowly if at all. Localization of the photolysis must 
be ascribed to low penetration by the exciting light. 

Existence of the photolysis has been ascribed! to disso- 
ciation of excited NOs;~ ions to NO.~ ions and O atoms. 
Under electron bombardment, the same process probably 
occurs; in addition, many NOs;~ ions will be stripped cf 
one outer electron, forming NOs radicals. This radical 
may later pick up an electron, with a probability of forming 
NO;- and O rather than NO;-; or it may decompose 
thermally to NOz and O. The NO:z may remain as such in 
the crystal, and hydrolyze on solution to HNO:+HNO;; 
or may pick up an electron to give nitrite ion. The small 
amounts of hydrogen gas and of alkali are ascribed to 
reaction with water of electrons which were trapped during 
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irradiation by barium ions, or at vacancies or imperfections 
in the crystal lattice. 

The authors are pleased to acknowledge the help and 
counsel of Professor Milton Burton, in charge of radiation 
chemistry at the Metallurgical Laboratory when this 
work was performed, and Dr. F. J. Safford, in charge of 
operating the Van de Graaff generator. 


* The work described herein was performed at the Metallurgical 
Laboratory of the University of Chicago under Contract No. W-7401- 
eng-37 with the Manhattan District. 

1 Narayanswamy, Trans, Farad. Soc. 31, 1411 (1935). 





Structure for Soap Micelles as Indicated by a 
Previously Unrecognized X-Ray 
Diffraction Band* 


R. W. Mattoon, R. S. STEARNS, AND W. D. HARKINS 
George Herbert Jones Chemical Laboratory, University of Chicago, 
Chicago 37, Illinois 
March 10, 1947 


EASUREMENTS on a previously unrecognized 
x-ray diffraction band By suggest that in an 
aqueous solution the micelle of a soap, or similar colloidal 
electrolyte, consists of only one double layer of soap mole- 
cules (Fig. 1). The subscript M designates that an asso- 
ciated spacing dy measures the approximate thickness of 
the micelle. Over a range of concentrations of a given soap 
dy is apparently constant and about twice the length of 
the particular molecule involved (Table I). 

McBain has considered that soap micelles exist in 
spherical highly ionic and lamellar weakly ionic forms; 
whereas, Hartley has considered them as spherical only. 
Evidence for the lamellar micelle as a stack of double 
layers of soap molecules, with “bound” water between the 
double layers, seemed to be given by the x-ray work of 
Hess, Kiessig, Philippoff, et al. (1937-1942). Similar experi- 
mental results have been discussed in two papers? from 
this laboratory in which extensive references to work on 
micelle structure are included. 

It is now found that the spacing dy’, associated with the 
new diffuse band, is essentially constant for any given 
soap, as it should. be if it represents the thickness of the 
double layer. The thickness increases from 36.5A for the 
8-carbon atom soap to 58.1A for the 16-carbon atom soap 
(Table I). This is an increase of 21.6A caused by the 
addition of 16 C-atoms to the thickness of the double 


TABLE I. Spacings in aqueous soap solutions calculated from the 
position of the peak of the micelle thickness band Byy. 








Bragg Modified Bragg 
spacing® 2 Xlength spacing» 
dy of molecule dy’ 
Potassium soap +1A +1A 





Caprylate, KCs 29.7 36.5 
Caprate, KCio 

Laurate, KCiz 

Myristate, KCis 

Palmitate, KCie 


Sodium lauryl! sulfates 








* dy =\/(2 sin@), valid for many diffracting layers. 

> dag’ =1.23X/(2 sin@), valid in the limiting case for two diffracting 
centers (first maximum of sinx/x). Note that dy <2l<dyj’ for the 
simple potassium soaps. 

* Kindly supplied by the Procter and Gamble Company. 
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Fic, 1. Cross section of a highly idealized micelle. A. Thickness ot 
micelle dyz’ ~50A for potassium myristate, distance between nearest 
hydrocarbon chains ds’ ~5.7A (Bragg spacing dg =dg’/1.23 ~4.6A, 
molecular area o¢ = (3)4dg/2/2 ~27A?2). B. Thickness of solubilized 
hydrocarbon layer Adyg’ ~ 13A for benzene. 


layer. The increase in double the length of the molecule 
is 82.54, or 20.3A, which is the same within the limit 
of error of the dy’ measurements. 

The short spacing ds, with a constant Bragg spacing of 
about 4.6A, is essentially independent both of the length 
of the soap molecule and of the absence or presence of oil 
solubilized in the micelle. This distance corresponds to a 
molecular area of 27A? in a plane perpendicular to the 
chains. 

The lowest soap concentration at which the micelle 
thickness band By has thus far been measured is 1.3 per- 
cent by weight of potassium palmitate. The band B,, 
related to the inter-micelle distance d;, begins to be notice- 
able above concentrations of about 9 percent potassium 
laurate and about 7 percent potassium myristate. The d,; 
spacing, which steadily decreases with removal of water, is 
always larger than dy’. The behavior of the three bands 
Bs, Bu, and B; with respect to position, intensity, and 
width, all as a function of the particular soap and its con- 
centration, seems to substantiate the following viewpoint. 
The micelle consists of one double layer of soap molecules 
from which the band By arises. As the concentration 
increases both the size and number of micelles probably 
increase. This means that the micelles are becoming closer 
together and a short-range order between micelles becomes 
more probable. This condition gives rise to the B; band. 

When a hydrocarbon is solubilized in a soap solution, the 
thickness dy’ of the micelle is increased an amount Ad,’ 
(Fig. 1). This gives additional evidence that dy’ represents 
the thickness of the micelle. For example, potassium 
laurate solutions saturated with ethyl benzene give a 
Ady’ of about 6A; potassium myristate solutions saturated 
with n-heptane give about 11A and saturated with benzene 
about 13A. These relations will be presented in detail later. 

In connection with this work the writers wish to thank 
Professor P. Debye for helpful advice. 


* This investigation was carried out under the sponsorship of the 
Reconstruction Finance Corporation, Office of Rubber Reserve, in 
connection with the Government's synthetic rubber program. 

1W. D. Harkins, R. W. Mattoon, and M. L. Corrin, J. Am. Chem. 
Soc. 68, 220 (1946). 

2W. D. Harkins, R. W. Mattoon, and M. L. Corrin, J. Colloid Sci. 
1, 105 (1946). 
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Sound Velocity and the Temperature Change 
of Molecular Association in Water 


ALFRED WEISSLER 
Naval Research Laboratory, Washington, D.C. 
March 10, 1947 





WO empirical criteria for molecular association in 

liquids, proposed recently by Lagemann,! are based 

on sound velocity. Unfortunately, one of these incorrectly 

predicts an increase in association with increase in tem- 
perature, and the other makes no prediction. 

Another quantity involving sound velocity appears 
better suited for providing information on temperature 
changes of molecular association. The molar sound velocity 
R= Mv!/d (where M is molecular weight, and v and d 
are sound velocity and density measured at the same tem- 
perature) was found by M. R. Rac? to be a temperature- 
independent constant for normal, non-associated liquids. 
But for water, and to a lesser extent such liquids as methyl 
alcohol and acetone, R increases considerably with tem- 
perature. It seems probable that R would remain constant 
even for these materials if the calculation utilized the 
effective M (which diminishes as association decreases) 
rather than the chemical-formula M. If this is true, in- 
creases in v!/d are a direct indication of decreases in extent 
of association. 

This is illustrated in Fig. 1, where the relative associa- 
tion (taken as 1.0000 at 0°C), found from the ratio of 
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Fic. 1. Relative association of liquids as a function of temperature. 


vo'/do to v:4/d:, is plotted as a function of temperature ¢. 
Water, which is highly associated, undergoes a decrease’ 
of about 7.2 percent in association from 0° to 100°; while 
toluene, which is not associated, shows nc appreciable 
change. Methyl alcohol occupies an intermediate position. 

Correlation with the results of other methods is difficult. 
Over a dozen techniques have been used*® to study the 


association of water, but the results disagree over a range | 


of hundreds of percent. One detailed investigation,® based 
on Raman spectra, indicates that the association of water 
decreases by about 13.2 percent from 0° to 98°. However, 
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molecular aggregates which are stable for only 10-4 second 
will show up in Raman scattering, but the aggregates 
must be stable for 10~* second in order to be detected by 
these sound velocity measurements. The latter method 
would therefore be expected to yield smaller values. 

It may be noted, further, that absolute values of associa- 
tion cannot be obtained from the ratio of observed to 
calculated molar sound velocities. This is true because the 
bond increments for calculating R already include the 
associating tendency due to, for example, the —OH group. 

1R, T. Lagemann, J. Chem. Phys. 12, 464 (1944). 

2M. R. Rao, Ind. J. Phys. 14, 109 (1 940) 

3 The values used, for sound velocity in Rites were those recently 
ae a at one megacycle in this Laboratory by Miss Patricia E. 
Morte 

4N. E. Dorsey, Properties of Ordinary Water Substance (Reinhold 
Publishing Corpcration, New York, 1940), p. 


5Q, Fruhwirth, Monats. f. Chemie 70, 157 (1037). 
6I. R. Rao, Proc. Roy. Soc. A145, 489 (1934). 





Non-Ionic Detergents as Association Colloids 
Giving Long X-Ray Spacings in 
Aqueous Solution 
JAMES W. MCBAIN AND SULLIVAN S. MARSDEN, JR. 


Department of Chemistry, Stanford University, California 
March 4, 1947 


OAPS and other ionic detergents in solution diffract 
x-rays so definitely as to demonstrate the presence of 
lamellar micelles. These lamellar micelles expand to take 
in layers of hydrocarbons when the latter are soiubilized 
in the detergent solutions. It has not been known whether 
such lamellar colloidal particles exist in solutions of non- 
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FIG, 2. 


ionic detergents. It is the purpose of this note to establish 
that non-ionic detergents also exhibit similar long spacings 
on x-ray examination in aqueous solution. 

The spacings observed are of the order of 40 to 150A 
and their value depends upon the concentration, as shown 
in the accompanying graphs. 

Para-isooctyl phenoxy nonacthoxy ethanol.—This is the 
condensation product of p-isooctyl phenol with ethylene 
oxide. Solutions of this material in water are all isotropic 
solutions exhibiting anomalous viscosity.! They produce 
x-ray diffraction patterns consisting of one or two orders 
of a very long spacing. There is also a broad halo due to 
side spacings with its maximum density at 4.5A. The long 
spacing of the solutions increases almost linearly with 
decrease in concentration, as shown in Fig. 1. This is 
similar to the now well-known behavior of ionic detergents 
in water. The liquid material itself does not give a pattern 
at room temperature but, when cooled to below its melting 
point, it gives a long spacing of 50A, which agrees well 
with the calculated value. 

Diethylene glycol monolaurate.—This surface active agent 
is supplied by Glyco Products Company, Inc. In aqueous 
systems at concentrations from about 35 percent to 85 per- 
cent, it forms two phases: one is isotropic and gives no 
long spacing; the other is liquid crystal and gives a long 
spacing which varies as the reciprocal of the concentration 
of the diethylene glycol monolaurate (see Fig. 2). This is 
decidedly in contrast with the behavior of ionic deter- 
gents. The material is not a pure chemical; the impurities 
are not completely known and so their effect on the long 
spacing has not been determined. 
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Glyceryl monolaurate.—This is supplied by Glyco Prod- 
ucts Company, Inc., as a white wax-like solid, and contains 
soap as a slight impurity. Its aqueous systems are also 
liquid crystalline and give long spacings which vary as the 
reciprocal of the concentration of the glyceryl monolaurate. 

Additional work is being done in this laboratory on the 
preparation and examination of pure non-ionic compounds 
of this type. This, together with a detailed description of 
the above work, will be reported elsewhere. 


1 Karl Baedeker, Kolloid. Zeits. 94, 161 (1941). 





Errata: Lattice Energies, Equilibrium Distances, 
Compressibilities, and Characteristic Fre- 
quencies of Alkali Halide Crystals 


[J. Chem. Phys. 5, 143 (1937)] 
Maurice L. HuGGins 
Kodak Research Laboratories, Rochester, New York 
March 7, 1947 


N columns 2 and 3 of Table I, M’ should be 12. 
In the fourth line of the second column on page 144, 


9 
the equation should be e= — qh? max 


In Eq. (3a), (1+ke) should be replaced by ke. 
In Eq. (3b), a in front of the bracket should be replaced 
by ke. 
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1 
The heading of column 7 of Table III should be (>>) ; 
B\OT/p 


The calculated compressibility values listed in columns 
20 and 21 of Table III are incorrect. The correct values, 
together with the experimental values from column 19, 
are given below. 
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(Eq. (2)) 
10-12 barye™! 


exp. 
(ref. 16) 
10-12 barye™! 
1.17 
3.41 


4.31 
6.01 


(Eq. (4)) 
Formula 10-12 barye71 
LiF 
LiCl 
LiBr 
Lil 


NaF 
NaCl 
NaBr 
Nal 


KCl 
KBr 
KI 


RbF 
RbCl 
RbBr 
RbI 


CsF 
CsCl 


CsBr 
CsI 
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